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This  report  describes  <33  investigations  in  the  field  of  seismic 
discrimination.  These  are  grouped  as  follows:  regional  seismic 
studies  (13  contributions),  miscellaneous  studies  (f>  contribu- 
tions), and  recent  developments  in  our  data  and  computer  systems 
(4  contributions). 
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This  is  the  twenty-eighth  Semiannual  Technical  Summary  report  describing  the  activities 
of  Lincoln  Laboratory,  M.I.T.,  in  the  field  of  seismic  discrimination.  This  report  covers  the 
period  1 October  1977  to  31  March  1978.  The  objective  of  the  Lincoln  Laboratory  program  is 
to  carry  out  fundamental  research  into  the  seismological  problems  associated  with  the  detection, 
location,  and  identification  of  earthquakes  and  nuclear  explosions,  and  with  the  estimation  of 
explosion  yields.  In  order  to  investigate  these  problems,  we  are  continuously  improving  our 
ability  to  manipulate  and  process  seismic  data  from  a global  network  of  high-quality  digital 
stations  and  arrays. 

During  the  present  six-month  period,  we  have  concentrated  mainly  upon  studying  the  nature 
of  seismograms  at  regional  distances  (<2000  km)  in  the  context  of  seismic  discrimination.  A 
variety  of  studies  have  investigated  the  utility  of  crustal  body-wave  phases  and  short-period 
surface  waves  in  the  detection,  location,  and  discrimination  of  events  at  these  distances. 

We  have  compared  locations  of  three  nuclear  explosions  in  Nevada,  using  crustal  phases 
only,  with  those  from  teleseismic  data  only  and  found  the  accuracy  of  the  epicenters  obtained  to 
be  roughly  equivalent,  but  that  crustal-phase  calibration  may  frequently  be  necessary.  The 
utility  of  differential  Pg-Pn  travel  times  as  a depth  discriminant  has  been  studied:  it  appears  to 
work  only  at  distances  of  less  than  500  km.  Crustal-phase  travel -time  data  from  NTS  explosions 
have  been  used  to  calibrate  the  Basm  and  Range  structural  province.  A study  of  film  chip  re- 
cordings cf  Soviet  explosions  has  demonstrated  that  although  Pn  and  Sn  are  readily  identifiable, 
the  other  crustal  phases  are  more  difficult  to  identify.  A comparison  of  broad -band  and  short- 
period  recordings  of  regional  events  at  the  Albuquerque  SRO  seems  to  indicate  some  advantages 
in  broad -band  recording  for  discrimination  purposes.  Two  sets  of  sonograms,  one  for  the 
Western  and  the  other  for  the  Eastern  U.S.,  show  the  viability  of  this  technique  in  phase  identi- 
fication. Polarization  filtering  > appears  to  enhance  phase  identification  at  regional  distances. 

The  Lg  phase  has  been  investigated  for  a variety  of  paths.  Wavenumber  spectra  at  NORSAR 
indicate  that  this  very  emergent  feature  of  regional  seismograms  is  mostly  incoherent  and  due 
to  scattering  from  more-coherent  higher-mode  Rayleigh  waves  of  short  duration.  The  propa- 
gation of  Lg  across  Western  Russia  has  been  studied  at  both  NORSAR  and  WWSSN  stations  in 
Scandinavia:  it  appears  to  be  transmitted  quite  efficiently  across  this  region,  with  no  anomalous 
attenuation  due  to  either  the  Ural  mountains  or  the  Baltic  Sea. 

Amplitude-distance  curves  for  various  crustal  phases  have  been  determined  for  Eastern 
North  America  and  applied  to  a larger  dataset,  including  Soviet  explosions,  for  magnitude  deter- 
mination. A magnitude  scale  based  on  Lg  amplitudes  appears  to  have  some  utility  as  a discrim- 
inant when  compared  with  body-wave  magnitudes. 

Several  studies  involving  teleseismic  data  also  have  been  carried  out.  A preliminary 
assessment  of  the  utility  of  ocean-bottom  seismometers  for  seismic  discrimination  purposes  has 
been  carried  out.  The  seafloor  appears  likely  to  provide  low-noise  sites  for  surveillance  of 
previously  inaccessible  regions,  and  recent  technical  developments  make  adequate  seismometer 
deployment  there  quite  feasible.  A separate  study  of  ocean-ridge  earthquakes  shows  that  their 
locations  can  be  considerably  improved  with  ray-tracing  procedures  in  a realistic  ocean-ridge 
model.  A method  for  the  extension  of  ellipticity  corrections  to  surface -reflected  phases  is 
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described.  A study  of  P-wave  velocity  in  the  Tonga  deep-earthquake  zone  shows  a correlation 
of  the  major  anomalous  velocity  region  with  a sudden  shear  strain  increase.  Analysis  of  Meshed 
SRC')  long-period  recordings  of  a double  explosion  indicates  that  events  separated  by  as  little  as 
20  sec  can  be  separated  in  the  surface -wave  train.  The  spectral  characteristics  of  events  in 
an  aftershock  sequence  have  been  studied  using  SRO  data:  seismic  scaling  effects  are  readily 
observable. 

We  continue  to  develop  our  data-handling  facilities.  A 300-Mbyte  disk  has  greatly  increased 
our  on-line  storage  capability.  System  software,  particularly  that  to  handle  Datacomputer  inter- 
action, has  been  substantially  improved.  The  seismic  display  package  has  now  been  implemented. 

M.  A.  Chinnery 


viii 


_ - • 


»— 


S K I S M I C DISCHIMINA'I  1 0 N 


I.  REGIONAL  SEISMIC  STUDIES 

A.  HIE  USE  OF  CRUSTAI  PHASES  IN  LOCATING  EXPLOSION  EPICENTERS 

We  relocated  the  Greeley.  Scroll,  and  Uulison  explosions  using  combinations  of  the  available 
data  set  of  arrival  times  in  an  attempt  to  evaluate  the  usefulness  of  crustal-phase  observations. 
Our  data  sets  for  Greeley  and  Scroll  consisted  of  all  the  arrival  times  reported  in  the  ISC  Bul- 
letin complemented  by  additional  crustal-phase  arrival  times  in  the  0 to  10°  range.  For  Hulison, 
we  had  only  the  ISC  data.  It  was  necessary  to  read  these  latter  arrival  times  from  the  LHSM  film 
because  there  were  no  S phases  listed  in  the  shot  reports.  A fortuitous  advantage  to  our  pro- 
cedure was  uniformity  of  observation,  since  all  the  events  were  read  by  one  analyst.  In  order  to 
use  all  the  data  for  each  event  in  a single  location,  we  assigned  the  errors  of  observation  given 
in  Tabl»  1-1  to  the  various  phases.  The  effect  of  this  was  to  reduce  the  dependence  of  the  epi- 
center on  arrival  times,  which  cannot  be  read  as  accurately  as  the  first  arrivals. 


TABLE  1-1 

ESTIMATED  ERRORS  OF  OBSERVATION 

Error 

Phase 

(sec) 

P 

0.5 

Pn 

1.0 

S 

1.0 

Sn 

1.0 

PKP(DF) 

1.0 

Pg 

3.0 

p* 

3.0 

Sg 

3.0 

S* 

3.0 

PP 

5.0 

PcP 

5.0 

The  ISC  epicenters  for  these  events  are  plotted  in  Fig.  I- 1 . Greeley  and  Scroll  were  825- 
and  20-kt  shots  detonated  at  NTS.  Rulison  was  a 40-kt  gas  stimulation  experiment  detonated  in 
Colorado.  Since  the  crustal  structures  vary  between  these  two  source  regions,  we  expect  to 
find  the  contribution  of  the  crustal-phase  arrival  times  sensitive  to  the  layer  parameters. 

The  epicenter  locations  were  done  using  the  Jeffreys- Bullen  travel-time  tables  for  tele- 
seismic  arrivals,  and  either  the  JB  tables  o^  our  own  Basin  and  Range  tables  (see  Sec.  C below) 
for  the  crustal-phase  arrivals.  In  each  case,  the  depth  was  constrained  to  that  given  in  the  shot 
report.  Jeffreys'  method  of  "uniform  reduction"  was  used  to  reduce  the  effect  of  erroneous  data 
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values.  IV  our  knowledge,  these  procedures  are  essentially  the  same  as  those  used  by  the  I'SGS 
and  ISC  Bulletins. 

The  results  of  a vari«  tv  of  location  experiments  are  given  in  Tables  1-2  to  1-4.  They  gen- 
erally tend  to  show  that  a large  amount  of  data  in  the  regional  distance  range  gives  epicenters 
as  accurate  as  those  derived  from  a large  amount  of  teleseismic  data.  On  the  other  hand,  add- 
ing a few  good  regional  stations  to  a teleseismic  data  set  usually  improves  the  resulting  epi- 
center both  in  accuracy  and  quality.  The  results  for  Hul.son  were  not  as  good  as  those  for 
Greeley  and  Scroll.  We  feel  that  this  is  due  to  our  having  to  rely  on  tabulated  crustal-phase 
arrival  times,  and  the  Colorado  plateau  crust  differing  from  either  our  Basin  and  Hange  struc- 
ture or  that  implied  by  the  JB  tables. 

The  following  general  conclusions  also  follow  from  our  data.  Of  the  six  possible  crustal 
phases,  only  Pn  is  reliable  enough  to  use  in  epicenter  locations.  The  others  are  useful  primar- 
ily for  determining  crustal  structure.  The  accuracy  of  crustal-phase  locations  is  directly  de- 
pendent on  how  well  the  crustal  velocity  structure  is  known.  In  unfamiliar  regions,  the  crust 
will  have  to  be  "calibrated"  in  order  to  use  near-in  stations  for  location  work.  Finally  several 
of  the  results  given  in  Tables  1-2  through  1-4  seem  to  show  high  quality  but  inaccurate  epicenters. 
These  appear  to  reflect  bias  in  the  travel -time  tables  which  should  be  eliminated  either  by  "cal- 
ibration" or  relative  event-location  techniques.  If  the  existence  of  a well-recorded  master 
event  in  a region  can  be  assumed,  then  substantial  improvements  to  our  results  may  be  possible. 

D.  W.  McGowan 
R.  K.  Needham 


B.  Pg  - Pn  THAVKI.  TIMES  FOB  NTS  EXPLOSIONS 


Figure  1-2  is  a plot  of  142  Pg  minus  Pn  travel  times  for  NTS  explosions  recorded  by  I HSM 
stations  mainly  in  the  Basin  and  Range  province.  The  two  lines  on  the  plot  are  the  theoretical 
zero  depth  Pg  - Pn  time  from  McCowan  and  Needham1  and  an  empirical  line  from  Ryall.2  The 
theoretical  line  is  based  on  the  following  plane  layered  Basin  and  Range  model: 


The  equations  are: 


Layer 

Thickness 

P-Velocity 

1 

15 

6.03 

2 

18 

6.73 

3 

an 

7.88 

t = t(Pg)  - t(Pn) 


J 


x2  + h2 


h tan  <p  i 


V 


3 


+ 


h 

Vj  cos  <p  j 


- A 


A = 


2 Dp  tan  4>p 


2Dj  tan^ 


2D. 


2D, 


+ Vp  COS  0p  + Vj  COS  (p  | 


= 6 sec 


where  si n<p^  = V^/V^  and  sin<|)p  = Vp/V^.  In  these  equations,  h is  the  source  depth  and  x is 
the  distance,  both  in  kilometers. 
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a. 


l'ln  partial  derivatives  evaluated  at  zfro  depth  are: 


in 

f'X 


h 0 


0 . 0 3 H 9 nec/km 


U\ 

Ph] 


h-0 


0.106  ui’c/km 


In  other  words,  this  (iiscriminant  is  almost  throe  times  more  sensitive  to  depth  than  it  is  to 
distance.  On  the  other  hand,  observational  errors  of  0.1  sec  limit  the  potential  depth  resolution 
to  ‘5  km-  Achieving  this  depth  resolution  w'mld  require  locations  to  be  known  within  1 km.  These 
estimates,  of  course,  assume  that  the  layer  parameters  are  constant  and  known.  Myall's  line^ 
is  based  on  NTS  explosion  and  Nevada  earthquake  data  out  to  a 300-km  distance.  Our  data  show 
substantial  numbers  of  explosions  above  the  line  in  the  250-  to  300-km  range.  Unfortunately, 
Myall's  line  falls  increasingly  away  from  the  data  population  beyond  300  km.  In  this  regard,  the 
theoretical  zero  depth  line  follows  the  data  population  but  also  illustrates  the  increasing  scatter 
with  distance. 

Our  conclusion  is  that  Pg  - Pn  travel  time  may  be  a uiscriminLnt  at  short  distances  but 
errors,  possib’y  structurally  dependent,  play  a dominant  role  beyond  500  km.  The  cause  of 
these  errors  should  be  studied,  so  that  they  can  be  eliminated  in  the  data  analysis. 

D.  W . McCowan 
M.  E.  Needham 


C.  CALIBR/  HON  OF  THE  BASIN  AND  RANGE  WITH  NTS 
CRUSTAL- PHASE  DATA 

We  reread  the  LRSM  short -period  film  to  enlarge  the  data  set  of  crustal- phase  travel  times. 
After  limiting  ourselves  to  20  NTS  shots  and  34  LRSM  stations  in  the  Basin  and  Range  (BR)  prov- 
ince, we  had  72,  20,  74,  39,  23,  and  32  Pn,  P*,  Pg.  Sn,  S*.  and  Sg  arrival  time::,  respectively. 
From  these,  we  determined  the  two-iayer  Basin  and  Mange  structure  given  in  Table  1-5.  Fig- 
ures 1-3  to  1-8  sho\  n?  corresponding  travel  times  fiom  our  model  plotted  along  with  the  data. 


BASIN  AND 

TABLE  1-5 

RANGE  AND  JEFFRE YS-BULLEN  MODEL  PARAMETERS 

BR  P 

JB  P 

BR  S 

JB  S 

Layer 

(km/sec) 

(km/sec) 

(km/Sec) 

(km/sec) 

■■ 

15.2 

15.0 

o.ul 

3.53 

3.37 

19.4 

18.0 

6.55 

3.73 

3.74 

00 

oO 

7.86 

4.43 

4.42 

it  appears  from  these  data  that  the  JB  crustal  model  is  an  excellent  approximation  to  the 
Basin  and  Range  except  in  the  first  layer.  Our  first-layer  velocities  are  8-  and  5-percent  higher 
for  P and  S,  respectively.  Furthermore,  the  Basin  and  Range  does  appear  to  be  well  modeled 
by  a two-layer  crust  as  far  as  crustal-phase  travel  times  are  concerned. 

D.  W.  McCowan 
R.  E.  Needham 
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L>.  CRUSTAL  l1  MASKS  FROM  SOVIET  PNEs  RECORDED 

AT  WWSSN  S TATIONS 

Wo  have  studied  WWSSN  short-period  recordings  in  the  distance  range  2°  to  16°  from  14 
PNl-iS  in  Western  Russia  in  an  attempt  to  determine  some  of  the  characteristics  of  propagation 
and  soui  ce  effects  in  this  area.  The  events  and  paths  used  are  shown  in  Fig.  1-9.  An  experi- 
enced analyst  picked  as  many  phases  as  possible,  without  recourse  to  travel-time  tables,  and 
the  resulting  suite  of  185  time  picks  from  35  seismograms  is  shown  in  Fig.  1-10.  We  attempted 
to  associate  seismic  phases  with  as  many  of  these  time  picks  as  possible,  using  as  a basis  a 
model  determined  by  Masse^  for  the  Canadian  Shield.  We  generated  travel-time  curves  for  the 
crustal  phases  Pn,  Sn,  P<,  S*,  Pg,  and  Sg  and  assigned  any  time  pick  falling  within  4 sec  (for  P) 
or  8 sec  (for  S)  of  the  gL  en  travel  time  to  the  particular  phase.  We  were  able  to  identify  35  Pn, 
32  Sn,  16  1’g,  13  Sg,  12  P*,  and  8 S*  in  this  manner.  A further  20  arrivals  were  identified  by 
the  analyst  as  the  phase  Lg  from  their  duration  and  character.  Thirty-nine  (21  percent)  of  the 
arrivals  could  not  be  identified  with  the  above  crustal  phases  and  were  probably  associated  with 
upper-mantle  refractions  with  apparent  surface  velocities  of  8.5  or  greater:  One  of  these 
(10.4  km/sec)  was  clearly  identifiable  and  corresponds  to  one  observed  by  Masstf  and  Alexander4 
in  the  same  region.  We  believe  that  our  results  show 

(1)  that  a shield  model  is  appropriate  for  the  paths  studied, 

(2)  that  Pn  and  Sn  can  almost  always  be  identified  anci  .r.easured, 

(3)  that  the  other  crustal  p’lases  are  more  difficult  to  identify,  and 

(4)  that  a traintd  analyst  can  pick  times  for  phases  prior  to  identifica- 
tion with  a false-alarm  rate  of  less  than  20  percent. 

For  every  time  pick  made,  an  amplitude  (peak-to-peak)  was  also  measured  except  where 
oif  scale,  and  reduced  to  ground  motion  at  the  station.  Foi  none  of  *he  events  were  there  suf- 
ficient stations  to  permit  construction  of  an  amplitude  distance  curve,  and  we  therefore  reduced 
toe  amplitudes  for  each  event  by  a scale  factor  related  to  ihe  ISC  through 

Reduced  Amplitude  = (Recorded  Amplitude)  X (5.0-m^) 

i.e.,  to  the  amplitude  from  an  event  of  mfa  = 5.0.  Amplitude-distance  curves  were  plotted  for 
each  phase:  the  results  wert  uniformly  hopeless  and  are  not  shown  here. 

For  comparison  with  the  other  studies  included  in  this  report,  we  have  determined  amplitude 
ratios  of  (Pn/Lg),  (Sn/Lg),  and  (Sn/Pn),  shown  in  Fig.  1-11.  Some  rather  dubious  identifications 
of  Lg  are  included  herej  in  some  caFes  only  an  upper  limit  to  its  amplitude  could  be  determined 
and  thus  only  a lower  limit  to  the  ratios  (Sn/Lg)  and  (Pn/Lg)  determined.  These  are  inoicated 
by  arrows  in  Fig.  I-H.  These  ratios  may  be  of  some  use  as  a discriminant,  but,  since  we  have 
studied  explosions  only,  they  are  included  here  purely  as  an  indication  of  the  nature  of  explosion 
seismograms  in  this  region.  Some*  of  these  ratios  have  been  included  in  a study  (elsewhere  in 
this  report)  of  explosion  and  earthquakes  in  the  Canadian  Shield,  or  similar  regions. 

R.G.  L'irth 
L.  C.  Lande 


E.  CRUSTAL  PHASES  RECORDED  ON  SRO  SHORT-PERIOD 
AND  BROAD- BAND  SYSTEMS 

Crustal  phase  Pn,  Pg,  and  Lg  recorded  on  the  broad-band  SRO  system  at  Albuquerque  are 
compared  with  the  corresponding  records  from  a short-period  vertical  SRO.  It  should  be  noted 
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TABLE  1-6 

EVENT  PARAMETERS  FOR  SRO  BROAD-P<\ND  RECORDED  EVENTS 

Date 

15  September  1977 

27  September  1977 

30  September  1977 

Origin 

14:36:30.1 

14:00:00.2 

10:19:19.6 

Latitude  (°N) 

37.033 

37.151 

40.51 

Longitude  (°W) 

116.043 

116.068 

110.28 

Depth 

n 

0 

5 

Mb 

4.5 

4.8 

5.0 

Location 

NTS 

NTS 

Utah 

Distance 

8.03 

8.07 

6.33 

Azimuth 

287.88 

288.67 

332.66 

Figure 

1-16 

1-17 

1-18 

that  the  records  came  from  separate  recording  systems  located  at  the  Albuquerque  Seismic 
Laboratory.  Record  pairs  from  two  NTS  shots  and  a Utah  earthquake  are  shown  in  Figs.  II-12(a) 
and  (b).  The  bottom  record  of  each  pair  is  from  the  short-period  SRO.  The  onset  of  the  Pn 
phase,  the  first  arrival,  from  first  explosions  failed  to  trigger  the  event  detector  for  the  short - 
period  SRO;  consequently,  this  record  is  not  complete.  Each  event  is  located  about  8°  from 
Albuquerque  with  the  great -circle  path  more  northwesterly  for  the  earthquake.  The  event  pa- 
rameters are  listed  in  Table  1-6. 

From  the  analog  records,  the  advantages  of  broad-band  recording,  if  any,  can  not  be  seen. 
However,  comparison  of  power  spectra  shows  that  longer-period  signals  that  are  preferentially 
excited  by  the  earthquake  are  suppressed  by  the  short-period  SRO  system.  This  can  be  seen  in 
Figs.  II-13(a)  and  (b)  by  comparing  the  Lg  spectra  in  the  bandpass  from  0.2  and  0.3  Hz.  The 
Lg  power  is  about  10  dB  above  the  noise  level  on  the  short-period  SRO,  anu  about  20  dB  above 
the  noise  level  on  the  broad-band  record.  This  disparity  is  consistent  with  the  12-dB/octave 
slope  in  the  amplitude  response  of  the  short-period  SRO  which  is  symmetrical  about  a peak  fre- 
quency of  1 Hz  (Ref.  5).  The  broad-band  Pg  signal  is  also  richer  in  longer-period  energy. 

All  power  spectra  were  computed  from  512  data  points  which  correspond  to  about  26  sec  of 
record  at  the  dip  itizing  rate  of  20  samples/sec.  Zeros  were  used  to  pad  out  the  time  windows 
of  signals  with  Durations  less  than  512  points.  Some  Pn  and  Pg  time  windows  (Figs.  1-12)  were 
significantly  shorter  than  26  sec,  with  the  result  that  ratios  of  signal-to-noise  power  inferred 
from  Figs.  1-13  and  1-14  are  underestimates  by  a factor  of  2 to  4.  A cosine  taper  was  used  to 
smocth  the  spectral  estimates. 

A comparison  of  the  power  spectra  from  the  NTS  shots  reveals  only  a marginal  increase  in 
longer-period  energy  in  the  broad-band  recordings.  Spectra  from  the  larger  of  these  shots  are 
shown  in  Figs.  I-14(a)  and  (b).  Spectra  from  a small  Colorado  earthquake  (not  shown  in  Figs.  1-12) 


reveah’d  significant  energy  above  the  Nyquist  frequency  of  10  Hz.  Aliased  energy  like  this  is 
• '\|)fcl ed  l e h»  a greater  problem  where  transmission  paths  are  more  efficient  such  as  those 
that  cross  shield  regions  or  stable  platforms. 

T.J.  Fitch 
R. M.  Sheppard 


F.  SONOGRAMS  OF  SRO  I1ROAD-RAND  RECORDED  EVENTS 

The  second  SRO  seismometer  located  at  the  Albuquerque,  New  Mexico  SRO  site  has  been 
modified  to  record  broad-band  data  directly  from  the  output  of  the  instrument.  These  signals 
are  the  raw  seismometer  outputs  and  are  not  filtered  in  any  manner.  Recording  of  the  data  is 
achieved  by  treating  the  signals  as  if  they  are  the  short-period  data  and  thereby  recording  the 
data  in  the  normal  position  on  tape  that  is  usually  occupied  by  the  sh  rt-period  data.  Both  the 
horizontal  and  th . vertical  outputs  are  recorded  at  20  H When  treated  in  this  fashioi  , the 
recording  and  analysis  of  the  data  can  be  done  with  most  of  the  existing  computer  programs. 

Several  well- recorded  events  have  been  analyzed  for  the  frequency  content  in  the  signal. 

Three  excellent  examples  of  these  events  will  be  discussed  here.  The  first  two  events  are 
underground  nuclear  explosions  at  the  Nevada  test  site,  and  the  third  event  is  an  earthquake 
located  in  Utah.  The  event  parameters  are  given  in  Table  1-6. 

Figure  1-15  shows  the  relative  distribution  of  the  event  locations  with  respect  to  ANMO. 
Figures  1-16,  1-17,  and  1-18  are  sonogram  plots  (frequency  vs  time)  of  the  explosions  and  the 
earthquake  described  in  Table  1-6.  In  these  figures,  the  vertical  axis  represents  frequency  and 
the  horizontal  axis  is  relative  time.  A sonogram  therefore  is  a representation  of  the  frequency 
content  of  a signal  vs  time.  To  achieve  these  figures,  the  event  signal  is  first  filtered  by  31 
constant-Q  filters,  then  the  31  outputs  are  integrated  over  1 -sec  intervals,  and,  finally,  the 
values  are  contoured  into  a sonogram.  The  black  areas  of  the  figures  represent  the  areas  of 
peak  signal,  while  the  hatched  areas  are  those  regions  which  are  6 dll  down  (1(  A amplitude  lower) 
from  the  peak  signal.  The  heavy  solid  black  lines  indicate  areas  which  are  12  dll  lower  than  the 
peak  areas.  Each  figure  displays  the  3 components,  Z,  N-S,  and  E-W,  of  the  ANMO  broad-band 
seismometer  (except  Fig.  1-16  which  did  not  record  the  N-S  component). 

The  areas  of  interest  in  Figs.  1-16  through  1-18  lie  in  the  frequency  distribution  of  the  sur- 
face waves.  The  two  explosions  exhibit  a more  confined  or  limited  frequency  distribution  than 
does  the  earthquake.  One  might  expect  that  the  explosion  would  contain  fewer  high  frequencies 
in  the  surface- wave  train  than  does  the  earthquake,  since  the  shots  are  2°  further  from  ANMO 
than  the  earthquake.  However,  there  seems  to  be  a wider  distribution  of  frequencies  for  the 
earthquake  than  one  would  expect  from  just  a 2®  increase  in  distance.  Tnis  wider  frequency 
distribution  is  also  evident  over  the  entire  waveform  of  the  earthquake.  It  is  evident  from  these 
figures  that  a discriminant  based  on  spectral  ratio  of  the  surface-wave  train  would  discriminate 
these  three  events.  However,  further  analysis  is  continuing  on  this  potential  discriminant  using 
broad-band  data. 

R.M.  Sheppard 

G.  SONOGRAMS  FOR  AN  EASTERN  U.S.  LRSM  PROFILE 
FROM  THE  NUCLEAR  EXPLOSION  SALMON 

We  analyzed  an  east  coast  profile  of  LRSM  3-component  short-period  recordings  of  the 
SALMON  explosion  by  the  sonogram  process.  The  explosion  and  station  locations  are  shown 
in  Fig.  1-19.  The  sonogram  output  is  a series  of  energy  contours  over  a l'requency-vs-time 
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TABLE  1-7 

SALMON  CRUSTAL-PHASE  DATA 


Station 

A 

(km) 

A 

(deg) 

Phase 

Travel 

Time 

Frequency 

(Sonogram) 

Frequency 
(Shot  Report) 

Amplitude 

(mp) 

EU-Al 

242 

2.18 

Pn 

36.5 

3.0 

4.0t 

370t 

P* 

39.5 

3.0 

- 

197 

Pg 

42.5 

3.0 

3.3t 

474t 

(SPZ) 

Sn 

67.5 

2.6 

- 

570 

(SPZ) 

Lg 

73 

2.0 

2.5t 

450 

BL-WV 

1058 

9.59 

Pn 

129.2 

2.6 

2.5?t 

5.1 2?  t 

Pg 

173.3 

2.1 

2.5t 

29. 2t 

(SPT) 

Sn 

257.1 

43 

(SPT) 

Lg 

300.7 

2.0t 

45 

BR-PA 

1375 

12.4 

p 

176.1 

2.8 

3.0t 

11. 5t 

(SPT) 

Lg 

388 

l.lt 

91. 2t 

(SPT) 

Sn 

301.5 

45 

DH-NY 

1795 

16.2 

p 

225.0 

2.1 

2.5f 

3.55t 

(SPT) 

Lg 

503 

I.:t 

66.4t 

LS-NH 

2121 

19. 1 

P 

265.5 

(SPT) 

Lg 

595 

l.ot 

49. 7t 

HN-ME 

2497 

22.5 

P 

298.3 

1.0 

l.lt 

17. 5t 

(SPZ) 

Lg 

705 

1.3t 

28 

(SPZ) 

Sn 

518 

15 

SV2QB 

3192 

28.7 

P 

358. 0t 

2.1 

2.0 

5.8t 

t Indicates 

values  taken  from  the  shot  repor 

t. 
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coordinate  system.  The  time  axis  is  equivalent  to  group  velocity  if  the  distance  to  the  event  is 
known.  In  our  implementation,  the  output  is  produced  by  applying  a suite  of  constant-Q,  narrow- 
band  filters.  Since  these  are  recursive  time-domain  operators,  there  is  a significant  group  de- 
lay between  arrivals  on  the  seismogram  and  energy  contours  on  the  plots. 

The  SALMON  explosion  was  detonated  in  Mississippi  and,  in  our  case,  observed  along  the 
whole  east  coast  from  2.18°  to  28.7°  distance.  It  is  thought  that  the  Eastern  U.S.  is  more  rep- 
resentative of  the  Eurasian  slueld  and  platform  areas  than  the  Western  U.S.  It  is,  at  least,  a 
high-Q  thick-crust  region  that  propagates  diagnostic  phases  like  Lg  in  contrast  to  the  Basin  and 
Kange  province  where  Lg  is  difticult  to  observe.  Unfortunately,  there  are  several  different 
geologic  environments  along  this  profile  with  different  crustal  structures.  The  explosion  and 
the  nearest  station,  EU-AL,  are  in  the  Gulf  Coast  plains.  The  intermediate  stations,  BL-WV, 
BH-PA,  and  Dll -NY,  are  on  the  eastern  plateau.  The  furthest  stations,  LS-NH  and  l IN -ME, 
are  on  the  New  England  plateau.  One  last  station,  SV2QB,  which  is  not  a part  of  the  profile  is 
on  the  Canadian  Shield.  Consequently,  we  expect  a substantial  scatter  in  observed  amplitudes 
and  periods  with  distance. 

As  noted  by  Alexander/*  the  obvious  advantage  to  the  sonogram  process  is  that  patterns  of 
energy  distribution  can  be  identified  which  would  be  lost  in  simple  amplitude  and  period  mea- 
surements. From  these,  empirical  discriminants  like  spectral  ratios  can  be  inferred  and  de- 
fined over  frequency-time  intervals.  Representative  sonograms  for  stat.  ns  ETJ-AL,  BR-PA, 
and  DH-NY  are  shown  in  Figs.  1-20  to  1-22,  respectively.  The  data  into1  1 for  station  EU-AL 
is  unique  because  it  includes  all  arrivals  from  the  first  P wave  through  the  surface-wave  train. 

The  results  obtained  by  means  of  these  energy  distribution  plots  are  listed  in  Table  1-7, 
along  with  some  taken  from  the  SALMON  shot  report.  We  were  able  to  identify  several  crustal 
body-wave  and  S-wave  phases  which  were  not  given  in  the  shot  report.  Sn-to-Lg  ratios  mea- 
sured off  the  seismograms  after  the  phases  were  identified  varied  from  1.25  at  EU-AL  to  0.54 
at  IfN-ME. 

Our  conclusion  is  that  the  sonogram  process  is  a useful  tool  for  identifying  seismic  phases 
from  +heir  energy  distributions  in  the  frequency- time  domain  and  that  it  should  be  applied  to 
data  from  more  regions  as  it  becomes  available. 

D.  W.  Me  Co  wan 

H.  POLARIZATION  FILTERING  OF  THE  FIRST  ARRIVALS 

FOR  THE  EASTERN  U.S.  LRSM  PROFILE  FROM  SALMON 

Polarization  filtering  has  Leen  used  to  aid  the  analysis  of  the  LRSM  profile  of  seismic  data 
from  the  U.S.  explosion  SALMON  described  above  in  Sec.  G.  Figures  1-23 (a)  and  I-24(a)  show 
the  3-component  data  from  EU-AL  to  HN-ME  ranging  from  A = 2.2°  tc  22.5°.  The  data  have  been 
rotated  toward  the  SALMON  epicenter,  and  phase-free  hi-pass  filtered  with  a 6-dB  cutoff  of 
0.8  Hz  to  i emove  microseisms.  The  times  have  been  reduced  in  the  figure  by  an  amount 
A/10  sec,  to  better  align  high-velocity  phases. 

A three-dimensional  polarization  scheme  was  used  to  enhance  planar  elliptical  motion,  e.g., 
isolated  surface  and  body  waves,  while  i ejecting  ellipsoidal  data  in  which  all  three  axes  are 
comparable.  The  least-squrxes  ellipsoid  of  motion  with  axes  l^>f2>f3  is  calculated  in  a 
moving  window  of  data,  and  the  two  largest  axes,  and  1 2,  are  projected  onto  the  data  com- 
ponents. This  amounts  to  truncating  the  data  expansion  to  the  two  most-significant  vector  com- 
ponents out  of  three,  and  results  in  a considerable  "cleaning  up"  of  3-component  data. 
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The*  projected  data  art'  further  scaled  by  a gain  factor  which  depends  on  the  instantaneous 
ellipticitv  t defined  by 

< vlf  ,/<!*  + *2]i  2 °'s  t 1 

where 


0 < <s'  1 

Since  f ^ < t,~  1 t,  t is  nearly  0 when  / ^ « i2-  This  case  occurs  when  the  data  are  polarized 
in  the  plane  of  an  ellipse  with  axes  and  1 2<  As  becomes  larger,  i approaches  1 in  the  case 

of  spheroidal  motion,  i.e.,  where  * 2 ~ * 3'  orc*er  to  rc.l°ct  spheroidal  motion,  in  the 

least-squares  sense  and  enhance  elliptical  motion,  a Butterworth-lype  ellipticity  filter  was  used 
as  the  gain  factor  given  by 

g - i [i  + (< A0)2n;,/Z  . 

This  gain  factor  is  a low-pass  filter  which  approaches  1 as  t — 0 and  decreases  On  dB/octave 
above  the  cutoff  value  <Q.  Parameters  n ard  are  arbitr  ary. 

In  Figs.  I -2 3(b)  and  I-24(b),  the  polarized  data  are  shown  displayed  under  the  original  data, 
^’or  this  processing,  the  cutoff  ellipticity  was  = 0.25,  the  filter  order  was  n = 3,  and  the 
mo\ing  time  window  for  the  ellipsoid  calcula*ions  was  1 sec  wide.  Comparing  these  figures, 
one  sees  that  the  polarization  filtering  sharpens  the  arrivals  and  extracts  a number  of  unknown 
phases,  compressional  and  shear. 

A detailed  analysis  of  all  the  enhanced  arrivals  is  not  possible  without  a clorely  spaced  net- 
work of  stations  as  used,  for  example,  in  the  Early  Pise  experiment.3 

Clearly,  polarization  methods  can  be  useful  for  near-regional  problems  in  discrimination, 
but  gcod-quality  3-component  data  must  be  obtained.  Detailed  analysis  of  complex  crustal  phases 
cannot  be  attempted  with  a sparse  network,  and  a dense  network  seldom  yields  unambiguous  geo- 
physical results. 

C.  W.  Frasier 

I.  TI'E  MECHANISM  OF  Lg  PROPAGATION 

Press^  lists  the  maximum  horizontal  Lg  to  maximum  vertical  Rg  amplitude  ratio  as  a pos- 
sible discrimination  parameter.  On  theoretical  grounds,  it  is  reasonable  to  expect  the  ratio  to 
be  significantly  greater  for  earthquakes  than  for  explosions  due  to  the  relative  distribution  of 
shear  and  compressional  energy  between  the  two  sources  and  the  suspected  nature  of  Lg.  In- 

g 

deed,  Bath  gives  the  measured  ratio  as  2 : 1 to  earthquakes  and  0.8  for  explosions.  These 
figures  have  not  been  established  for  statistical  quantities  of  events,  and  earthquakes  that  violate 
the  discriminant  can  be  found.  The  problem  may  lie  in  our  poor  understanding  of  the  nature  of 
Lg  - its  generation,  its  propagation,  and  its  observed  structure.  To  better  understand  Lg,  its 
spatial  and  temporal  behavior  over  NORSAR  subarrays  has  been  investigated. 

Briefly,  Lg  is  observed  over  continental  source-receiver  paths  as  a short-period  large - 
amplitude  wave  with  an  average  source-receiver  velocity  of  about  3.5  km/sec.  Most  of  its 
motion  is  transverse,  though  an  appreciable  vertical  component  is  present.  The  phase  is  ex- 
tremely emergent,  often  building  up  over  several  minutes  before  peaking,  l igure  1-25  shows  the 
the  seismogram  recorded  on  the  NORSAR  short -period  vertical  (SPV)  instrument  at  location  01  BO 
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for  the  event  listed  in  the  Pl)K  as  4 October  1971  at  10;00;02.0  at  61.6"N  and  47.1  K (Western 
Russia!  of  magnitude  9.1  m^.  The  event  is  presumed  to  be  an  underground  explosion.  It  mea- 
sures 17.4"  from  NORSAR  with  a great -circle  back  azimuth  of  65".  'The  bottom  scale  shows 
source-receiver  group  velocities.  At  a velocity  of  3.5  km/sec  (marked  by  Lg),  there  is  a local 
maximum  following  a continual  buildup  of  energy  beginning  approximately  3 min.  before.  The 
data  lie  in  the  band  0.5  to  2 IIz.  Figure  1-26  snows  the  records  at  each  element  of  the  IB  sub- 
array. While  there  is  coherence  in  the  envelope,  small-scale  coherence  is  low. 

q 

Thi>  wavenumber  spectra  of  several  time  sections  of  these  records  were  examined  to  de- 
termine the  portion  of  coherent  energy  and  the  direction  and  phase  velocity  of  its  propagation. 

Th  > intervals  studied  were:  the  first  arrival.  1 min.  after  the  first  arrival,  2 min.  prior  to 
tiie  maximum  Lg  time,  at  the  maximum  Lg  time,  immediately  after  the  maximum  Lg  time,  and 
1 min.  past  the  Lg  time.  Of  these  intervals,  only  the  firnt  arrival  and  the  section  immediately 
after  the  maximum  Lg  time  (but  not  including  it)  showed  any  appreciable  wavenumber  coherence 
in  the  source  direction.  The  wavenumber  spectra  for  these  two  periods  are  shown  in  Figs.  1-27 
and  1-28  Accordingly,  the  vertical  component  of  the  Lg  phase  for  this  event  as  recorded  at 
NOHSAR  may  be  describee  as  follows;  a steady  buildup  of  essentially  spatially  incoherent  energy 
peaking  at  a group  velocity  of  3.5  km/sec,  followed  by  a coherent  wave  train  of  time  duration  less 
than  25  sec  arriving  from  *he  direction  of  the  source  with  phase  velocity  4.1  km/sec,  followed  by 
diminishing  incoherent  motion. 

The  above  observations  are  consistent  with  the  interpretation  of  the  vertical  part  of  Lg 
(Ref.  9)  as  second-  and/or  higher-mode  Rayleigh  waves,  with  the  bulk  of  their  energy  confined 
to  a part  of  the  upper  crust  (but  not  at  the  surface),  and  the  emergent  and  diminishing  parts, 
including  the  maximum  trace  amplitude,  being  scattered  waves.  The  consequence  of  such  be- 
havior is  that  the  amplitude  in  the  scattered  portion  of  the  phase  is  not  only  a function  of  dize 
of  the  propagating  wave  but  also  of  the  scale  of  inhomogeneity  in  the  crust  near  the  receiver. 

The  amplitude  of  the  propagating  part  of  the  phase  would  be  tne  integr  ated  effect  of  spreading 
and  scattering  over  the  entire  path.  Under  thes^  circumstances,  the  usefulness  of  Lg  for  re- 
gional detection,  location,  and  discrimination  purposes  would  be  enhanced  by  employing  small 
arrays  and  spatial  processing  that  could  separate  path,  source,  and  receiver  characteristics 

T.  E.  Landers 

J.  Lg  MEASUREMENTS  AT  NORSAR  FROM  SOVIET  EXPLOSIONS 

AND  TWO  EARTHQUAKES 

We  have  investigated  the  available  short -period  data  from  NORSAR  in  an  attempt  to  cali- 
brate Lg  propagation  across  Western  Russia,  and  to  determine  if  topographic  features  such  as 
the  Urals,  the  Parent  Seas,  and  the  Baltic  Sea  are  barriers  to  this  propagation.  These  data 
complement  another  study  by  R.G.  North  (Sec.  D above)  which  used  W'W'SSN  stations,  and  the 
results  obtained  here  are  similar.  Table  1-8  is  a list  of  the  usable  events  which  had  long 
enough  data  windows  after  the  first  P arrivals  to  include  detectable  Lg  phases.  All  events 
except  those  with  nonzero  depth  are  explosions. 

Figure  1-29  shows  the  ISC  epicenter  locations,  m^  values,  and  dates  of  the  events  on  which 
we  measured  phases.  Seismograms  of  the  events  were  searched  for  emergent  Lg  arrivals  with 
group  velocities  between  about  3.5  to  3.6  km/ sec.  The  envelope  peaks  were  measured  where 
the  energy  was  clearly  above  the  noise  before  the  P arrival,  and  isolated  by  obvious  changes 
in  noise  and  frequency  character  before  and  after  the  Lg  phase.  Although  the  amplitudes  arc 
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NORSAR 

TABLE  1- 

EVENTS  FOR  WHICH  Lg  TIME 
EXPLOSIONS  HAVE 

8 

WINDOWS  WERE 
ZERO  DEPTH 

AVAILABLE; 

Event 

Date 

Location 

(°N)  (°E) 

mb 

■ 

Delta 

Lg 

Seen 

Lg  Amplitude 

i 

12/22/71 

47. 87  48. 22 

6.0 

0 

24.95 

Barely 

388 

2 

10/04/71 

61.61  47.22 

4.6 

0 

17.35 

Yes 

105 

3 

04/16/72 

47.79  16.04 

4.7 

18 

13.^0 

Yes 

664 

4 

06/17/72 

48.40  14.52 

4.6 

20 

12.63 

Yes 

64 

5 

07/09/72 

49.78  35.40 

4.8 

0 

17.73 

Yes 

108 

6 

09/30/73 

51.66  54.54 

5.2 

0 

25.6 

Good 

195 

7 

10/26/73 

40.  84  48. 29 

4.4 

41 

30.27 

No 

- 

8 

08/14/74 

68. 94  75. 83 

5.4 

0 

27. 45 

Yes 

80 

9 

08/29/74 

67.23  62.  10 

5.0 

0 

22.74 

Yes 

78 

10 

09/13/74 

49.  76  78.  03 

5.2 

0 

38.  11 

Maybe 

<25  (noise) 

Earthquakes  3 


TABLE  1-9 

NORSAR  TiME-DOMAIN  MEASUREMENTS 


P (biggest  arrival  in  first  20  sec) 

— L 9 

Distance  Amplitude  Period  Amplitude 

(deg)  (mp)  (sec)  (mp) 


Lg/P 

Amplitude 

Ratio 


0.9  0.32 

0.41 


probably  accurate  to  within  a factor  of  2 across  the  NOHSaH  array,  it  must  be  emphasized  that 
arrival  times  are  very  suspect  and  cannot  be  considered  accurate  for  location  purposes  to  within 
±30  sec  in  most  cases. 

Two  PNEs  in  Western  Russia  and  two  Austrian  earthquakes  are  closest  to  NORSAR.  The 
initial  P waves  of  the  two  PNEs  are  much  richer  in  high  frequencies  than  the  P waves  of  the 
Austrian  earthquakes.  Although  the  Austrian  earthquakes  are  4°  to  5°  closer  to  NORSAR  than 
the  PNEs,  the  structural  paths  are  more  complex  and  it  is  uncertain  how  much  of  the  high- 
frequency -explosion  energy  is  simply  passed  by  the  high-Q  Baltic  shield  and  Russian  platform. 
The  Lg  phases  for  all  four  events  are  very  emergent  and  complex  in  character  with  little  co- 
herence from  sensor-to-sensor.  This  visual  lack  of  coherence  occurs  even  across  the  sensors 
of  individual  subarrays  at  NORSAR.  Using  high -resolution  wavenumber  analysis,  however, 
Landers  has  separated  the  scattered  from  the  direct  waves  in  tne  Lg  waveform  (see  Sec.  I above). 

Time -domain  measurements,  picked  on  single  NORSAR  sensors  for  these  four  events  are 
summarized  in  Table  1-9.  We  note  that  the  two  earthquakes  have  inconsistent  P and  Lg  am- 
plitudes even  though  the  mb  and  locations  are  nearly  identical.  This  is  probably  due  to  different 
source  radiations  along  the  takeoff  direction  to  NORSAR.  The  PNE  amplitudes  seem  more  con- 
sistent even  though  their  locations  are  widely  separated.  Lg/P  amplitude  ratios  are  1.9  to  4.5 
for  the  earthquake  pair,  whereas  the  PNE  ratios  are  very  consistent  — 0.32  and  0.41. 

In  Fig.  1-30  the  spectra  of  a 20-sec  P window  and  a similar  length  Lg  window  are  displayed 
for  the  Austrian  earthquake  (event  3),  which  has  the  largest  Lg/P  amplitude  ratio.  Figure  1-31 
shows  similar  spectra  for  a PNE  (event  5).  The  high-frequency  character  of  the  explosion 
P wave  and  the  low  amplitude  of  the  Lg  spectrum  verify  the  time -domain  observations. 

Measurements  of  Lg  amplitude  for  these  four  events  ancl  the  remainder  of  the  PNEs  in 
Fig.  1-29  are  plotted  vs  distance  from  NORSAR  in  Fig.  1-32.  The  amplitudes  have  been  normal- 
ized to  mb  6.0.  The  dashed  trend  line  is  Nuttli's  theoretical  curve ^ for  Lg  decay  with  distance 
for  earthquakes  in  the  Eastern  U.S.,  which  besides  geometrical  spreading  includes  attenuation 
with  distance  given  by  exp[-0.07A|,  where  A is  measured  in  degrees.  This  curve  has  been 
shifted  vertically  to  pass  through  the  level  of  scaled  Lg  amplitudes  for  the  PNEs.  The  solid 
curve  above  the  data  is  Nuttli's  empirical  formula  for  Lg  amplitude  vs  distance  from  an  m^6.0 
earthquake.  This  comes  from  the  relation  between  short-period,  surface-wave  amplitude  and 
mb  given  by 

mb  = 1.66(log  A)  + log  (A/T) 

where  A is  measured  in  millimicrons  on  a vertical  short-period  instrument. 

From  Fig.  1-32  we  see  two  interesting  trends  of  the  Lg  data.  First  of  all,  there  is  no  clear 
effect  on  Lg  transmission  of  any  topographic  features  in  Fig.  1-29  such  as  the  Ural  Mountains  or 
the  Barent  Sea.  The  variation  of  Lg  amplitude  with  distance  from  NORSAR  seems  to  fit  Nuttli's 
theoretical  trend  for  the  Eastern  U.S.  Secondly,  the  general  level  of  Lg  phases  for  PNEs  when 
scaled  to  mb  6.0  have  a trend  at  least  a magnitude  unit  below  that  expected  for  earthquakes  in 
the  Eastern  U.S.  This  is  a short-period  observation  analogous  to  the  mb:Mg  criterion  which 
works  effectively  for  large  events  at  teleseismic  distances.  Clearly,  many  earthquakes  from 
Russia  are  needed  to  investigate  the  use  of  Lg  as  a discriminant.  The  data  from  the  two  coin- 
cident Austrian  earthquakes  are  quite  ambiguous  and  remind  us  that  short-period  surface  waves 
are  probably  as  sensitive  to  source  radiation  effects  as  are  short-period  body  waves. 

C.  W.  Frasier 
R. M.  Sheppard 


K.  '’HOF  U'.  ATION  OF  THE  Lg  PHASE  ACROSS  Till-:  URALS  AND  THE  BALTIC 


Ur  have  attempted  to  determine  to  what  extent,  *f  at  all,  natural  features  such  as  the  Urals 
at  I tin*  Ualtie  affect  Eg  amplitudes.  Because  of  the  vertical  absence  of  earthquakes  providing 
source-receiver  paths  contained  within  Western  Russia  and  the  apparently  similar  Baltic  shield 
regions,  we  have  studied  short-period  recordings  at  WWSSN  stations  ip  Scandinavia  from  three 
nuclear  explosions.  These  are  situated  to  the  west  of  the  Urals  (I’NE  shot  on  4 October  1971), 
south  of  the  Urals  (Western  Kazakh,  23  December  1969),  and  one  east  of  the  Urals  (Eastern 
Kazakh,  30  November  1969).  These  events,  their  recording  stations,  and  the  source-receiver 
paths  art  shown  in  Fig.  1-33.  Epicentral  distances  range  from  11  to  40°.  Paths  to  UME,  COP, 
and  KON  cross  the  Baltic;  those  to  KEV  and  NUR  do  not.  Lg  was  readily  observed  from  these 
large  explosions  even  up  te  40  distance;  its  amplitude  was  generally  quite  small  on  the  vertical 
component,  and  we  have  measured  amplitudes  on  the  horizontal-component  recordings  only.  The 
largest  peak-to-peak  amplitude  within  30  sec  of  the  time  corresponding  to  the  3.5-km/sec  group 
velocity  was  measured.  The  Lg  amplitude  was  in  all  cases  not  much  smaller  than  that  of  F;  its 
dominant  period  (1  to  2 sec)  is  unfortunately  close  to  that  of  the  prevailing  noise. 

Figure  1-34  shows  horizontal-component  Lg  amplitudes  as  a function  of  distance  A,  corrected 
for  peak  magnification  and  adjusted  to  6.0  through 

Amplitude  (6.0)  Amplitude  * 10^’° 

Note  that  these  results  seem  to  imply  that,  for  an  explosion  of  = 3.0,  the  amplitude  of  Lg  will 
be  at  least  of  amplitude  1 mp  at  distances  up  to  30".  Short-period  detection  levels  as  small  as 
1 mp  have  been  claimed:  thus,  if  Lg  is  of  any  use  in  a discrimination  context  it  should  be  mea- 
surable for  extremely  small  explosions  across  similar  paths. 

Amplitudes  measured  on  the  north -south  component  were  always  larger  than  those  measured 
on  the  east-west  component.  The  paths  shown  on  Fig.  1-33  run  nearly  in  an  east-west  direction, 
and  thus  the  maximum  amplitude  of  Lg  appears  to  be  in  the  transverse  component  of  motion, 
i.e.,  of  Love-wave  type. 

For  all  the  events  studied,  amplitudes  at  stations  to  the  west  of  the  Baltic  are  very  close  to 
those  observed  to  the  east.  The  Baltic,  an  extremely  shallow  sea  within  only  a few  isolated  pock- 
ets of  depths  greater  than  200  m,  thus  appears  to  have  no  observable  effect  upon  the  propagation 
of  Lg. 

Also  shown  on  Fig.  1-34  is  a curve  determined  empirically  by  Nuttli10  for  rhe  diminution  of 
Lg  amplitudes  with  distance  in  Eastern  North  America,  given  by 

A(a)  = k(A  sin  a)"  1//2  e'0,07A  (I-D 

4 

where  A is  distance  in  degrees,  and  k is  a constant,  here  set  to  10  , This  appears  to  be  a good 
fit  to  the  data:  note  that  the  Lg  amplitudes  measured  for  the  Eastern  Kazakh  event  are  quite  con- 
sistent with  Nuttli's  relation.  This  may  well  indicate  that  the  Urals  have  little  effect  upon  Lg 
amplitudes. 

Nuttli1®  has  also  derived  an  empirical  relation  between  my  and  Lg  amplitudes,  given  by 

mb  = 1.66  logA  + log  (A/T)  (1-2) 

wherp  A is  the  peak-to-peak  amplitude  in  millimicrons  of  Lg,  and  T is  its  dominant  period. 
Assuming  T = 1 sec.  Fig.  1-35  shows  m^  calculated  from  Lg  as  a function  of  distance.  The 


16 





avenge  m.  indicated  is  about  5.0,  or  one  m.  less  than  that  measured  from  teleseismic  P-wave 
n b 

observations.  In  addition,  Nnttli  used  Lg  amplitudes  measured  from  the  vertical-component 
seismogram  only,  which,  as  we  mentioned  earlier,  were  considerably  less  than  those  observed 
on  the  horizontal  components  from  which  we  measured  A.  It  is  thus  possible  that  Lg  amplitudes 
may  be  useful  as  a discriminant,  if  Eq.  (1-2),  determined  for  Eastern  North  America,  is  ap- 
plicable to  Western  Russia. 

R.G.  North 


L.  Pn.  Sn,  AND  Lg  AMPLITUDE-DISTANCE  RELATIONSHIPS 

FOR  EASTERN  NORTH  AMERICA 

Amplitudes  of  the  crustal  phases  Pn,  Sn,  and  Lg  for  the  Eastern  North  American  events 

listed  in  Table  I- 10  were  read  from  short-period,  vertical-component  seismograms  of  the 

Canadian  seismograph  network.  As  an  aid  in  identification  of  the  Pn  phase.  Massy's  P-wave 

velocity  model  was  used  to  calculate  theoretical  travel  times.  Calculations  of  theoretical  Sn 

wave  velocities  assumed  - 3 V 

P s 

Figures  1-36  through  1-38  present  amplitude-distance  plots  normalized  to  mfa  = 3.6  for  the 
three  phases.  Though  all  show  clear  trends  to  amplitude  decay  with  distance,  Pn  data  are  the 
most  widely  scattered  for  distances  of  less  than  10°.  Small  events  6 and  8 with  Pn  magni- 
tudes of  3.9  introduce  this  scatter.  Examination  of  Fig.  1-37  reveals  that  Sn  amplitudes  for  these 
same  events  are  mo’  F'nsely  clustered  about  the  line  of  average  decay  than  the  Pn  amplitudes 
of  Fig.  1-36.  The  ^ ..ter  consistency  of  Sn  readings  is  attributable  to  larger  amplitudes  (1.5  to 
2 times  Pn  amplitude)  at  the  same  distances.  These  observations  suggest  that  the  scatter  in 


TABLE  1-10 

EPICCNTRAL  DATA  FOR  EVENTS  USED  IN  THIS  STUDY 


Event 

Date 

Time 

Latitude 

(°N) 

Longitude 

(°W) 

Magnitude 

mb 

1 , Baffin  Island 

4/24/73 

03:55:04 

71.7 

71.2 

4.7 

4.2 

2.  Quebec- Maine 

6/15/73 

01 : 04:05 

45.3 

71.2 

5.0 

4.8 

3.  Blue  Mountain  Lake 

7/15/73 

08:20:31 

44.0 

74.4 

3.4 

- 

Blue  Mountain  Lake 

7/15/73 

10:32:38 

44.0 

74.4 

3.2 

3.4 

Mt.  Wright,  P.Q $ 

1/1 7/7  A 

15:58:02 

52.7 

67.6 

3.5 

- 

5.  Minnesota 

7/09/75 

14:54:15 

45.6 

96.0 

- 

4.6 

6.  Maniwaki,  P.Q. 

7/]2/75 

12:37:15 

46.5 

76.2 

- 

4.3 

7.  Western  Quebec* 

6/21/75 

21:42:51 

48.1 

78.0 

3.6 

- 

8.  Southwestern  Quebec 

1V19/75 

15:25:11 

47.0 

78.8 

3.8 

- 

9.  Atlantic  Ocean 

\0/06/75 

22:21:41 

49.6 

57.0 

- 

5.2 

t Mine  blast 

J Rockburst  in  Malarctic  Mine 
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Pn  amplitudes  (for  small  events)  may  be  due  to  comparatively  low  signal-to-noise  ratios,  and 
that  magnitudes  derived  from  Pn  amplitudes  for  such  events  will  be  less  useful  in  discrimination 
than  magnitudes  derived  from  the  Sn  and  Lg  phases. 

The  line  of  average  decay  for  Pn  amplitudes,  shown  in  Fig.  1-36,  was  found  by  Everndm.1* 

A method  similar  to  his  yielded  an  average  value  for  Sn  decay  with  distance.  The  least-squares 
fit  to  a straight  line  was  found  for  log^  (A/T)  against  logjQ  (distance)  for  each  event  with  suf- 
ficient data.  Six  events  combine  to  give  an  average  value  of  —1.80,  represented  as  a straight 
line  in  Fig.  1-37.  The  same  procedure  gives  an  average  value  for  Lg  amplitude  decay  of  —1.65. 
This  number,  the  average  of  eight  events,  agrees  closely  with  Nuttli's*^  value  of  —1.66.  Fig- 
ure 1-38  shows  a straight  line  with  a slope  of  —1.66  as  well  as  a curved  line  representing  the 
theoretical  attenuation  of  dispersed  surface  waves.  The  equation  for  this  line  is: 

A - K(A_1//3)  (sinA)‘1/2  e‘vA 

where  A is  the  Airy  phase  amplitude,  K is  a constant,  A is  the  distance  in  degrees,  and  y is 
the  coefficient  of  anelastic  attenuation. 

Table  I- 11  lists  the  events  used,  the  values  of  Sn  and  Lg  amplitude  decay  with  distance,  and 
the  standard  deviations  of  the  data  points  from  the  least-squares  fit  for  each  event.  In  Sec.P 
below,  these  amplitudes  will  be  converted  to  magnitudes  for  a comparison  between  earthquake 
and  explosion  values. 

M.  W.  Shields 
T.  J.  Fitch 
R.  E.  Needham 


TABLE  1-11 

AMPLITUDE  FALLOFF  WITH 

DISTANCE 

Event 

Sn  Slope 

Standard 

Deviation 

Lg  Slope 

Standard 

Deviation 

#15/73 

-1.74 

0.20 

-1.31 

0.17 

7/15/73a 

- 

- 

-1.58 

0.07 

7/1  V73b 

- 

- 

-1.52 

0.15 

1 1/1  #74 

-1.94 

0.15 

-1.93 

0.12 

7/09/75 

-1.98 

0.24 

-1.84 

0.26 

7/1  #75 

-1.73 

0.23 

-1.48 

0.15 

#21/75 

-1.61 

0.24 

-2.01 

0.17 

1 #19/75 

-1.77 

0.14 

-1.55 

0.26 

Average 

-1.80 

0.20 

-1.65 

0.17 
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M.  COMPARISON  OF  CRUSTAL- PHASE  MAGNITUDES  FOR  EXPLOSIONS 
AND  EARTHQUAKES 

Magnitude  relationships  of  the  usual  form 

m a + b log  A f logA/T  (1-3) 

are  used  to  test  for  systematic  differences  between  explosions  and  shallow  earthouakes  recorded 
in  the  near-regional  distance  range  from  5*  to  20°.  In  the  range  the  crustal  phases  Pn,  Sn,  and 
Lg  are  prominent  of  short-period  seismograms.  Magnitude  formulas  pertinent  to  Pn  and  I.g  were 
taken  from  Evcrnden1 1 and  Nuttli,*  ° respectively.  The  one  forSn  was  inferred  from  the  amplitude- 
distance  relationship  determined  empirically  by  Shields,  Fitch,  and  Needham  in  Sec.  L above. 

The  study  yielded  -1.80  for  the  coefficient  of  the  log  A term  in  Eq.  (1-3).  A value  of  3.79  for  the 
normalizing  constant  a in  Eq.  (1-3)  was  determined  by  equating  teleseismic  m^  to  mgn  for  five 
earthquakes  for  which  we  have  measured  Sn  amplitudes  and  periods.  Figure  1-39  represents  a 
test  of  the  Sn  amplitude-distance  relationship.  The  test  is  successful  in  that  the  magnitudes  at 
individual  stations  scatter  more  or  less  uniformly  about  the  mean  Sn  magnitude  for  the  event 
regardless  of  the  distance  between  source  and  receiver.  The  three  magnitude  relationships  are 
as  follows: 

mgn  = 3.79  + 1.80  log  A + logA/T 

mLg  = 3*30  + 1.66  log  A + log  A/T  (Nuttli10) 

mPn  + 2*0°  log  A 1 logA/T  (Evernden**) 


where  A is  in  degrees,  and  A/T  is  in  microns /second. 

Magnitudes  were  determined  for  one  mine  blast  in  Eastern  Canada,  four  PNEs  in  Western 
Russia,  the  SALMON  shot  in  Mississippi,  and  ten  earthquakes  located  in  Eastern  North  America. 
Amplitudes  and  periods  were  read  from  analog  short-period  records  from  LRSM,  WWSSN,  or 
the  Canadian  network.  All  magnitudes  refer  to  maximum  amplitudes  of  the  vertical  component 
of  motion.  The  various  magnitude  estimates  are  listed  in  Table  1-12.  Figure  1-40  chows  that 
the  ratio  of  Sn-to-Lg  amplitude  for  explosions  and  earthquakes  from  this  1imited  and  possibly 
heterogeneous  data  set  defines  distinct  populations  in  the  distance  range  from  600  to  1300  km. 
Removing  the  distance  dependence  by  converting  amplitudes  to  magnitudes  does  not  change  the 
basic  result.  The  m£g  from  explosions  is  about  0.3  of  a magnitude  unit  less  than  m^'s  for 
earthquakes  (Table  1-12).  Sn  magnitudes  for  explosions  and  earthquakes  are  in  close  agreement 
with  P-phase  magnitudes,  either  mpn  or  m^.  Consequently,  the  magnitude  ratio  mLg/mSn  will 
be  distinct  for  these  explosions  and  earthquakes  as  shown  in  Fig.  1-41. 

From  these  results  a reasonable  working  hypothesis  would  be  that  expl?sions,  either  be- 
cause of  their  shallow  depth  or  the  nature  of  the  source,  are  comparatively  poor  exciters  of  Lg. 
Apparently.  Sn  excitation  from  explosions  and  earthquakes  is  comparable  to  P-wave  excitation 
from  these  sources.  For  this  to  be  so,  the  Sn  from  explosions  must  be  primarily  a P-to-S  re- 
flection at  the  free  surface  near  the  source.  Furthermore,  the  efficient  propagation  of  Sn  or 
Pn,  and  their  respective  apparent  velocities,  implies  that  these  phases  are  trapped  in  the  litho- 
sphere, as  would  be  the  case  for  wide-angle  reflections.  This  model  for  Sn  and  Pn  propagation 
is  also  consistent  with  the  amplitude-distance  plots  for  these  phases  shown  in  Sec.  L above.  Pn 
amplitudes  show  no  pronounced  distance  dependence,  whereas  Sn  amplitudes  do.  Apparently,  the 
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TABLE  1-12 

MAGNITUDE  ESTIMATES  FOR  EARTHQUAKES  AND  EXPLOSIONS 


Earthquakes 


Event 

mb  (Sn) 

mb  (Pn) 

mb  (ISC) 

Mn 

4/24/73 

4.5  (7)a 

4.7  (3) 

4.2 

4.8  (7) 

V15/73 

5.0  (6) 

4.9  (5) 

4.8 

4.9  (6) 

7/15/73 

3.2  (3) 

3.2  (1) 

- 

3.4  (ISC) 

7/15/73 

- 

- 

3.4 

3.2 

1 1/30/73 

4.3  (3) 

4.5  (1) 

- 

4.5  (6) 

7/09/75 

4.6  (6) 

4.4  (6) 

4.6 

4.8  (7) 

7/\2/75 

4.3  (8) 

3.7  (6) 

4.3 

4.2  (8) 

6/21/75 

3.7  (4) 

3.8  (2) 

- 

3.6  (5) 

10/0^75 

5.2  (7) 

5.3  (7) 

5.2 

- 

1 Vi  9/75 

3.8  (4) 

3.9  (3) 

- 

3.8  (4) 

Explosions 

lQ/2i/64b 

4.4  (5) 

4.3  (23) 

4.4C 

4.0  (10) 

l/lV74d 

3.9  (7) 

- 

- 

3.5  (7 

7/1Q/71® 

4.7  (3) 

4.8  (3) 

5.2 

4.2  (3) 

9/19/71 

4.1  (2) 

4.5  (4) 

4.5 

3.7  (3) 

1C/04A1 

4.4  (3) 

5.0  (3) 

4.6 

3.9  (3) 

9/04/72 

4.4  (3) 

4.3  (4) 

4.6 

4.2  (3) 

a.  Numbers  in  parentheses  are  those  of  amplitude  readings. 

b.  SALMON  nuclear  test. 

12 

c . From  Jordan  et  a]_. 

d.  Western  Quebec  mine  blast 

. 

e.  All  following  explosions  were  located  in  Western  Russia. 

J 
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signal -to-noisr  ratio  for  tin'  Pn  signals  is  significantly  lower  than  that  for  the  Sn  signals.  This 
disparity  can  L>e  explained  in  terms  of  the  wide-angle  reflection  hypothesis.  Pn  energy,  having 
the  higher  phast  velocity,  tends  to  leak  energy  out  of  the  lithosphere,  whereas  Sn  energy  tends 
to  be  trapped.  Of  course,  scattering  and  anelastic  attenuation  will  contribute  to  amplitude  dim- 
inution of  both  phases.  This  speculation  should  be  tested  with  more  observations  and  by  com- 
parisons of  real  and  synthetic  seismograms. 

T.  J.  Fitch 
M.  W.  Shields 
II.  E.  Needham 
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Fig.I-1.  Event  epicenters  for  Western  U.S.  explosions. 
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Fig.  1-6.  Sn  travel  times  for  NTS 
explosions  and  best-fitting  model 
(39  observations). 
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Fig.  1-7.  S*  travel  times  for  NTS 
explosions  and  best-fitting  model 
(23  observations). 


Fig.  1-8.  Sg  travel  times  for  NTS 
explosions  and  best-fitting  model 
(32  observations). 
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Fig.I-11.  Amplitude  ratios  of  (Pn/Sn),  (Pn/Lg),  and  (Sn/Lg) 
measured  from  seismograms  as  a function  of  distance. 
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Fig.  1-12 (a).  Broad-band  and  short-period  SRO  recordings 
of  Pn  and  Pg. 
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Fig.I-12(b).  Broad-band  and  short-period  SRO  recordings  of  Lg. 
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Fig.I-13(a).  Utah  earthquake,  Fig.I-13(b).  Utah  earthquake, 

30  September  1977,  vertical  30  September  1977,  SP-Z  SRO 

broad-band  SHO  at  AIjQ.  at  AL/Q. 
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Fig.I-14(a).  NTS,  27  September  1977,  Fig.  I-14(b).  NTS,  27  September  1977, 

SP-Z  SRO  at  AIjQ.  vertical  broad -bend  SRO  at  ALQ. 
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Fig.  1-15.  Distribution  of  events. 
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Fig.  1-19.  Hast  coast  LRSM  stations  for  SALMON 
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Fig.  t-<£1.  S \ ; MON  sonograms,  station  MK-PA, 
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Fig.  1-23.  (a)  Rotated  3-component  data  from  SALMON  recorded  at  LRSM 

stations  LU-AL,  BL-WV,  and  HR -PA.  (b)  Polarized  data  as  described  in 
text  with  cutoff  ellipticity  t = 0.25  and  filter  order  n = 3. 
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tig.  1-29.  Map  showing  epicenters  of  presumed  explosions  in  Western 
Russia  and  three  earthquakes  recorded  at  NORSAR.  These  events  have 
sufficiently  long  recorded  codas  to  include  Lg  arrivals  on  Sp  vertical 
sensors. 
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Fig.  1-30.  Spectra  of  20-sec  windows  of  P,  Lg, 
and  noise  before  P of  event  3 — Austrian  earth- 
quake of  16  April  1972.  Note  that  signal/noise 
ratio  of  P phase  is  nearly  20  dBfrom  1 to  4 Hz. 
mb  = 4.7,  A - 13.5°. 
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Fig. 1-31.  Spectra  of  20-sec  win- 
dows of  P,  Lg,  and  noise  before 
P of  event  5 - PNE  of  9 July  1972. 
Note  higher  frequency  character  of 
P wave  and  lower  le^el  of  Lg  spec- 
trum compared  with  Pig.  1-30. 

•4.8,  A 17. 7°. 
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Fig.  1-32.  NOHSAH  measurements  of  Lg 
amplitudes  normalized  to  m^  6.0.  This 
was  done  for  each  event  by  scaling  Lg 

amplitude  by  factor  10^"mb.  Dashed 
curve  is  Nultli's  theoretical  trend  of  Lg 
amplitude  vs  distance  for  earthquakes  in 
Eastern  E.S.,  shifted  hero  to  pass  through 
data.  Solid  curve  is  Nuttli's  empirical 
variation  of  Lg  with  distance  for  m^  6.0 
earthquakes  based  on  a short-period, 
surface -wave  magnitude  scale. 
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WEST  LONGITUDE  (dtg) 

Fig.  1 — 33.  Source-receiver  paths  used  in  this  study. 
Urals  are  indicated  by  shaded  region. 
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Fig.  1-34.  Lg  amplitudes,  separated 
into  measurements  from  short-period 
north -south  and  east -west  component 
seismograms.  Also  shown  is  Nuttli's10 
relation  for  amplitude  dimension  of  Lg, 
with  k 104. 
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Fig.  1-35.  mb  calculated  from  Lg 
using  Nuttli's10  relation. 
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Fig.  1-36.  Pn  amplitude  vs  distance. 
Numbers  correspond  to  event  numbers 
in  Table  1-10.  Slope  of  line  of  average 
decay  is  —2.0. 


Fig.  1-38.  Lg  amplitude  vs  distance. 
Slope  of  straight  line  is  —1.66.  Co- 
efficient of  anelastic  attenuation  for 
theoretical  curve  is  y = 0.07  deg“*. 
Numbers  correspond  to  event  num- 
bers in  Table  I -10. 


Fig.  1-37.  Sn  amplitude  vs  distance. 
Numbers  correspond  to  event  numbers 
in  Table  1-10.  Slope  of  line  is  —1.80. 


iV-  * 

'•  !*\  M 
? 2 

4 ^ >. 
AV. 

*•  \\  & 
8.  W 


u 0« 


o t 4 e a io  it  u w >•  io  i* 


delta  (dag) 

Fig.  1-39.  Sn  magnitude  vs  distance. 
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II.  MISCELLANEOUS  STUDIES 


A.  OTEAN-HOTTOM  SEISMOMETERS  FOR  NUCLEAR  MONITORING 

AND  RESEARCH:  A REASSESSMENT 

Ocean-bottom  seismic  stations  may  play  an  essential  role  in  a comprehensive  nuclear  test 
ban  treaty  (CTBT)  which  allows  only  a limited  number  of  land  stations  to  be  operated  by  a moni- 
toring nation  within  a host  nation.  This  stems  from  the  likelinood  that  ocean-bottom  stations, 
located  in  international  waters,  will  not  be  subject  to  regulation.  They  would,  in  effect,  be  free 
stations  which  could  be  used  to  monitor  all  coastal  and  island  areas.  Significantly,  the  world's 
most  tectonically  active  areas  (see  Fig.  II- 1),  where  evasion  schemes  are  likely  to  be  attempted, 
are  near  the  deep  ocean.*  Clearly,  internal  land  stations,  whose  number  and  location  are  likely 
to  be  restricted  by  a host  nation,  are  not  necessary  to  monitor  these  high-seismicity  areas.  In 
any  event,  the  available  land  stations  should  not  be  wasted  here. 

The  use  of  ocean-bottom  stations  for  regional  surveillance  is  predicted,  of  course,  on  the 
ability  of  such  stations  to  detect  and  discriminate  small  explosions  from  earthquakes.  Unfor- 
tunately, adequate  information  is  not  available  to  assess  the  relative  performance  of  ocean- 

bottom  stations  in  a near-in,  or  regional,  context.  Initial  experiments  in  the  early  1960's  sug- 

2-5 

gested  that  the  ocean  bottom  may  be  a high  SNR  environment.  However,  further  studies  spon- 
sored by  ARPA  in  the  mid-60's  indicated  that  ^an-bottom  seismographs  (OBSs)  would  not  be 

6 **  8 

as  effective  as  land  stations  for  teleseismic  surv  llance.  This  conclusion  was  based  largely 
on  the  observation  that  the  lower-frequency  (<5-Hz)  background  noise  level,  measured  at  several 
ocean  sites  by  short-period  instruments,  was  much  higher  than  that  at  land  sites.  Accordingly, 
the  concept  of  using  ocean-bottom  seismographs  for  teleseismic  monitoring  was  abandoned  in  the 
late  1960's.  Today,  such  a judgment  may  no  longer  be  valid  in  view  of  the  substantial  advances 
in  seismic  instrumentation  and  oceanographic  methods  for  deploying  and  recovering  seafloor 
devices. 

9 

We  have  recently  completed  a preliminary  review  of  OBS  and  marine  technology.  From 
this  work,  it  appears  that  modern  digital  seismometers  rigidly  emplaced  beneath  the  seafloor 
by  deep-sea  drilling  ships,  manned  subrnersibles,  or  remote-controlled  manipulator  ships  could 
provide  low-noise  ocean-oottom  stations  useful  for  both  teleseismic  and  regional  monitoring.  In 
fact,  by  incorporating  modern  acoustic  data  telemetry,  it  should  be  possible  to  have  essentially 
real-time  monitoring  comparable  to  that  from  land  stations.  The  more  important  scientific  ad- 
vantages of  OBS  stations,  particularly  sub-seafloor  installations,  are  summarized  in  Secs.  1 
and  2 below.  Figure  II -2  shows  the  primary  seismic  ray  paths  on  ocean  bottom. 

1.  Seafloor  OBSs  vs  Land  Seismograph  Stations 

(a)  Increased  Signal  Amplitude:-  Although  most  previous  OBS  experiments  have  prob- 
ably not  measured  the  true  ground  motion  on  the  seafloor  because  of  poor  coupling  and  high  back- 
ground noise,  it  is  clear  that  signal  amplitudes  from  earthquakes  observed  on  the  ocean  floor  are 
generally  higher  than  those  seen  by  a land  station.  For  example,  the  ARPA-sponsored  Texas 

Instruments  (T.I. ) field  tests  off  the  Aleutian  Islands  showed  that  OBS-calculated  m,  values  av- 

7 b 

eraged  0.2  unit  greater  than  land -station  calculated  values.  This  was  thought  to  result  from  the 
fact  that  rays  traveled  a slightly  shorter  path  to  the  ocean-bottom  stations  and,  more  important, 
they  did  not  propagate  through  a low-Q  continental  crust.  Also,  more-recent  work  suggests  that 
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shallow  occar.c  paths  will  show  e.\tremelv  low  attenuation  compared  with  continental  paths.*  ^ 
Oceanic  (■)  .a lues  topically  exceed  oOOO  at  20  Hz. 

(b)  Lower  Noise:  i)eep*ocean  sites  are  relatively  isolated  from  the  sources  of  back- 
ground noise  likely  to  alfect  seismic  stations,  namely:  cultural  noise,  breaking  surf  on  coast- 
lines, storm  microseisins,  and  local  sea  surface  waves  of  tidal  current  effects.  In  practice, 
this  appears  to  be  the  case.2,5,7  The  dominant  noise  source  is  believed  to  be  surf  microseisms 
propagating  out  from  the  coastline  as  llayleigh  waves  in  the  water  mass  and  along  the  water- 
seafloor  interface  [Fig.  II-2(a)]. 


TABLE  ll-l 

COMPARISON  OF  SRO  STATION  NOISE  BACKGROUND 
WITH  OBS  OBSERVATIONS 
[Mean  rms  Noise  Amplitude  (mp)] 

vertical  Component 

Operational 

Dates 

Stat  i ory'lnvesti  gators 

Short  Period  (f  ) 

0 

Long  Period  (fQ) 

SROs+ 

Albuquerque,  N.M. 

0.38  (2.86  Hz) 

9.56  (25  sec) 

1977- 

Guam  Island 

40.25  (2  86  Hz) 

11.25  ( 25  sec) 

1977- 

New  Zealand 

28.92  (2.86  Hz) 

45.92  (25  sec) 

1977- 

OBSs 

Asada  and  Shimamura^ 

^50  (2.7  Hz) 

- 

1972- 

Francis  et  al  J2 

Francis  et  al  J3 

~25  to  50  (8.0  Hz) 

- 

1972- 

Texas  Instruments  Inc.^ 

*100  to  300(1  Hz) 

- 

1965-1968 

Sutton  et  al  p 

- 

5,000  (15  sec) 

1965-1972? 

Bradner  et  ol 

~800  (1  Hz) 

« 10, 000  (5  to  10  sec) 

1964 

Ewing  and  Ewing^ 

*1  (2  Hz) 

- 

1959-1961 

t Source:  Texas  Instruments  Quarterly  Report  J 1 

Table  1 1 - 1 lists  representative  noise  levels  observed  on  land  and  on  the  seafloor.  Note 
that  only  the  early  Lamont  workers2  reported  short -period  noise  amplitude  levels  comparable 
to  those  observed  on  land  instruments  (i.e.,  1-mp  p-p  in  the  2-  to  10-Hz  band).  Unfortunately, 
their  low-noise  observations  were  not  substantiated  in  the  more-extensive  field  tests  carried 
out  by  T.I.  in  the  mid-1960's.^  Noise  amplitudes  were  typically  two-order3-of-magnitude 
greater  than  those  reported  earlier  (i.e.,  100  to  300  mp.  at  1 Hz).  However,  it  must  be  empha- 
sized that  these  apparent  high-noise  observations  may  not  have  been  a true  measure  of  solid 
earth  motions  after  considering  the  methods  used  to  make  background  noise  observations.  In 
T.I.'s  system  and,  for  that  matter,  with  most  systems:  (1)  the  seismometers  simply  rested  on 
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surficial,  unconsolidated  se;  floor  sediments  and  ooze  whose  seismo-acoustic  properties  were 
not  much  different  than  the  overlying  water  mass  (Fig.  II-2(b):  OBS  position  1],  and  (2)  these 
devices  used  tall  vertical  frames  which  protruded  into  water  mass  to  house  their  seismometers. 
These  facts  suggest  that  most  OBSs  have  probably  recorded  ocean-water  motion  as  well  as 
solid-earth  ground  motion.  Significantly,  those  seafloor  OBSs  with  seismometers10,12,13  well- 
coupled  to  the  solid  earth  and  isolated  from  water-induced  motions  on  the  housing  frame  have 
shown  short-period  noise  levels  approaching  the  1-  to  10-m|ji  levels  observed  at  the  better  la  id 
SRO  stations  today  .'Table  II - 1 ).  No  studies  of  long-period  seafloor  noise  other  than  the  pre- 
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liminary  Lamont  and  Scripps  work,  ' have  been  reported  to  date. 

2.  Sub-seafloor  vs  Seafloor  OBSs 

Additional  advantages  can  be  gained  for  OBSs  by  implanting  the  seismometer  beneath  the 
seafloor.  Some  of  these  are  enumerated  below. 

fa)  Better  Coupling  to  the  Solid  Earth:-  The  most  obvious  method  for  improving  the 
coupling  of  ocean-bottom  seismometers  to  the  solid  earth  and  isolating  them  from  disturbing 
noise  in  the  ocean  water  and  surficial  sediments  is  to  rigidly  mount  the  seismometers  beneath 
the  seafloor  (Fig.  II - 2 (b ) : CBS  position  2].  Competent  sedimentary  roc’’3  usually  lie  within  a 
few  tens  of  meters  beneath  the  seafloor.  In  fact,  the  hard  igneous  rocks  of  the  high-Q  oceanic 
crust  are  only  covered  by  a few  hundred  meters  of  sediment  in  most  areas.  A borehole-type 
OBS  installation  on  or  within  the  oceanic  basement  (positions  3 or  4)  should  prove  vastly 
superior  to  the  conventional  free-fall  devices  previously  deployed.  The  signal  amplitude  can 
be  maximized  by  adjusting  the  seismometer  case  density  to  improve  its  acoustic  impedance 
match  with  the  surrounding  borehole  materials. 

(b)  Lower  Noise:-  A buried  seismometer  will  be  well-isolated  from  local  noise  gen- 
erated in  t’  e ocean  water.  The  air-water  and  seafloor-water  interfaces  serve  as  efficient 
reflectors  which  effectively  trap  c»ny  propagating  waves  in  the  water  [Fig.  II-2(c ) ).  In  fact, 
recent  tests  in  shallow  water  showed  more  than  a factor-of-10  reduction  in  local  ocean- 

14 

generated  noise  on  a vertical  component  seismometer  buried  only  2 m beneath  the  seafloor. 
Also,  in  a deep-water  test  conducted  aboard  the  GLOMAR  CHALLENGER  drilling  ship,  very 
low  noise  levels  («120  mp  at  10  Hz)  were  observed  beneath  the  seafloor.15 

The  distant  coherent  noise  propagating  in  the  water  mass  and  surficial  sediments  should 
also  be  markedly  attenuated  at  depths  of  only  a fev  undred  meters  since  the  shear-wave 
velocity  of  the  surficial  unconsolidated  sediments  is  only  about  0.1  to  0.4  km /sec  (Ref.  16). 

This  noise,  which  is  generated  beneath  the  breaking  waves  on  beach  surfaces,  will  remain 
primarily  in  the  waveguide  formed  between  the  reflecting  seafloor  and  underlying  harder-rock 
interfaces  (Fig.  I(-2(a)].  Seismometers  buried  300  to  600  m beneath  the  seafloor  would  be 
effectively  shielded  from  this  dominant  source  of  seaflocr  noise. 

(c)  Reduced  Signal  Contamination;-  Seismic  signals  received  at  a seafloor  OBS 
travel  through  the  ocean  water  as  well  as  through  the  solid  earth  beneath  the  station  (position  1 
in  Fig.  II -2 ).  Rays  reflected  from  the  local  air-sea  surface  interface,  as  well  as  shear  waves 
propagating  along  the  soft  sediment  waveguide,  contaminate  the  direct  seismic  rrival  phases. 
These  arrivals  not  only  introduce  apparent  complexity  in  the  wave -train  cod'*  Lut  they  also 
tend  to  reduce  signal  amplitude  because  of  inter' 'rjnee  effects.  By  locating  the  seismometer 
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in  a borehole  beneath  the  seafloor  surface,  seismic  waves  entering  the  overlying  water  mass 
and  soft  sediments  w:l'  be  effectively  trapped  much  like  the  noise  initially  generated  in  the 
ocean  (positioi  3 4). 

(d)  More  - form  Crust  and  Mantle  Structure:  The  crust  and  mantle  structure 
beneath  the  ocear.  \ as  . ns  is  now  known  to  be  much  simpler  than  that  beneath  the  continents. 

The  seafloor  spreading  hypothesis,  generally  accepted  by  ocean  scientists  to  account  f->r  the 
formation  of  the  seafloor,  predicts  that  nearly  horizontal  rock  layers  underlie  most  jf  the 
deep  ocean.  Marine  geophysical  measurements  support  this  hypothesis. 

Simple  layering  implies  that  large-aperture  arrays  could  be  deployed  to  further  improve 
the  SNR  by  beamforming  or  velocity  filtering.  In  fact,  many  of  the  signal-processing  problems 
caused  by  the  near-field  complexity  of  earth  structure  at  land  arrays  (e.g.,  LASA,  NORSAR) 
will  not  be  encountered.  By  sharply  reducing  signal-generated  coherent  noise,  a closer  ap- 
proach to  the  ideal  nTn"  SNR  enhancement  might  be  realized.  The  widespread  uniformity  of 
earth  structure  beneath  the  ocean  also  implies  that  much  larger  arrays  than  those  practical 
for  land  installation  could  be  built. 

J.  D.  Phillips 
D.  W.  McCowan 


B.  RELOCATION  OF  OCEAN-RIDGE  EARTHQUAKES 

Detailed  knowledge  of  the  location  of  ocean-ridge  earthquakes  relative  to  each  other  and  to 
the  ridge's  islands,  bathymetric  features,  and  magnetic/gravity  fields  would  be  useful  for  geo- 
physical studies  as  well  as  nuclear  monitoring  activities.  Unfortunately,  conventional  methods 

for  determining  epicenter  locations  on  ridges  may  not  be  adequate  to  attain  the  required  preci- 

17  18 

sion  and  accuracy  in  that  relatively  simple  Jeffreys -Bullen  or  Herrin  spherical  earth 

continental-type  velocity  models  are  used  to  calculate  travel  times  and  depths. 

The  computation  of  the  location  of  earthquakes  which  are  recorded  by  many  observatories 
distributed  in  an  azimuthally  uniform  pattern  around  the  epicenter  is,  of  course,  insensitive 
to  the  errors  in  any  radially  symmetric  earth  velocity  model  since  a travel-time  error  along 
one  ray-path  azimuth  is  compensated  by  errors  in  the  opposite  sense  along  reciprocal  azimuth 
ray  paths.  However,  few  oceanic-ridge  earthquakes  are  ever  teleseismically  recorded  by 
such  an  ideal  station  distribution  due  to  their  generally  small  magnitude  (m^  ~ 4.0  to  5.0). 
Typically,  only  those  stations  along  a nearby  continental  margin  record  the  smaller  events 
(mb  < 4.5).  This  fact  coupled  with  the  known  higher-velocity,  2 -dimensional  structure  beneath 
ridges  at  shallow  depths  (<200  km),  compared  to  continents,  means  station  azimuth-dependent 
errors  are  ikely.  Epicenter  location  and  depth  errors  can  be  expected  to  be  biased  along  the 
azimuth  toward  the  greatest  number  of  recording  stations. 

This  biasing  effect  of  epicenter  positions  was  recently  demonstrated  for  the  central  mid- 

Atlantic  ridge  by  comparing  the  epicenter  distribution  with  the  known  tectonic  axis  along  median 

19 

valley  and  transform  faults.  A markedly  skewed  arrangement  showing  wide  scatter  was  ob- 
served. The  epicenters  were  generally  located  to  the  north  and  west  of  the  tectonic  axes.  This 
distribution  was  attributed  to  the  "slowness"  of  the  JB  model  since  most  of  the  recording  sta- 
tions were  located  to  the  northwest  in  North  America. 

A more-realistic  earth-velocity  model  has  recently  been  devised  to  account  for  both  the 
higher  velocity  and  two-dimensional  nature  of  ocean-ridge  structure.  For  the  ridge  crest  re- 
gion dewn  to  a depth  of  78  km,  the  velocity  structure  proposed  by  Steinmetz  et  al.^  was  used. 
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This  structural  model  which  includes  shallow -depth  velocity  inversions  is  based  on  seismic 

refraction  results  from  the  mid-Atlantic  ridge  between  45°  and  40  °N  latitudes  and  petrologic 

considerations.  For  the  deep-ridge  structure  and  the  continental  regions,  the  Herrin  velocity 
1 8 

model  was  used.  The  ridge  axis  water  depth  was  taken  as  3 km.  Hay  tracing  was  then  done 
to  calculate  travel  times  and  epicentral  distances  for  earthquakes  occurring  on  the  ridge  axis 
at  depths  of  3,  33,  96,  and  200  km  (Refs.  21  and  22).  Figure  11-3  shows  comparisons  of  the 
differences  in  travel  times  between  the  "fast"  oceanic-ridge  model  proposed  here  and  the 
"slow"  ,)H  model.  Note  that  the  largest  differences  are  for  the  surface  depth  event.  This 
might  be  expected  since  the  ridge  model  shows  very  high  velocity  layers  (Vp  * 7.7  to  8.0  km/sec) 
within  5 km  of  the  surface.  In  fact,  true  mantle-like  layers  (Vp  > 8 12)  are  only  79  km  below 
the  surface.  The  JB  model  puts  these  layers  at  33  and  160  km,  respectively. 

To  test  the  efficacy  of  the  model  to  better  locate  earthquakes,  ten  events  which  occurred 
along  well-charted  portions  of  the  north  mid-Atlantic  ridge  during  an  arbitrary  period  (Mareh- 
September  1%4)  were  analyzed  (Fig.  11-4).  The  median  valley  and  transform  fault  axes  here 
are  known  with  an  accuracy  of  ±2  km.  Table  II -2  summarizes  the  relocation  parameters  of 
epicenters  shown  in  Fig.  II -4  as  well  as  for  other  nearby  events. 

Careful  examination  of  Table  II -2  reveals  that,  for  those  location  solutions  which  were 
depth  convergent,  the  new  depths  of  the  three  larger-magnitude  events  (>4.8)  became  more 
shallow  while  the  epicenter  moved  toward  the  ridge  axis  ( Fig.  II -4 ).  In  contrast,  the  smaller- 
magnitude  events  showed  increased  depths;  however,  their  epicenters  also  moved  closer  to 
the  ridge  axes.  Since  most  oceanic -ridge  earthquakes  are  small  and  are  believed  to  occur 
above  the  base  of  the  lithosphere  ( — 110  km  thick),  a test  was  made  to  determine  the  maximum 
epicenter  relocation  effect  by  depth  constraining  all  solutions  to  an  arbitrary  depth  of  113  km. 

The  open  circles  in  Fig.  II -4  show  these  relocations.  Note  that  in  all  cases,  save  one  (ISC 
event  829).  the  epicenters  moved  toward  the  ridge  axis  and  in  several  cases  they  moved  ex- 
actly onto  the  axis.'  Although  these  results  suggest  the  new  velocity  model  better  approximates 
the  true  velocity  structure  beneath  ridges  than  does  the  JR  model,  further  statistical  evaluation 
of  the  general  applicability  of  the  mode1  is  required.  Relocation  of  a large  number  of  events 
over  a more  extensive  portion  of  the  ridge  crest  would  be  appropriate. 

J.  D.  Phillips  L.  C.  Lande 

R.  E.  Needham  D.  W.  McCowan 

R.  S.  Sheppard 

C.  A NOTE  ON  ELLIPTICITY  CORRECTIONS  FOR  SURFACE-REFLECTED  PHASES 

23 

Dziewonski  and  Gilbert  using  addition  theorem  for  I,egendre  functions  presented  an  effi- 
cient and  accurate  method  of  calculation  of  ellipticity  corrections. 

They  proposed  to  evaluate  the  ellipticity  correction  according  to  the  formula: 

2 

<5t  = V rm(A)  P2  m(cosrQ)  cos  m4  (II-l ) 

m=0 

where  r is  the  colatitude  of  the  source,  i is  the  azimuth  from  the  source  to  the  receiver,  A is 
o 

the  epicentral  distance,  and 
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RELOCATED  EPICENTERS  FOR  OCEAN-RIDGE  EARTHQUAKES 
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rm,A)  P_1  ,j3  n?  am(0)d0  ^ ^m(0)  f(T,Z 

* Y K^noi2  -p2)1/2))k 

k 


where 

V0)--Tr,2.m<cos0> 

and  the  remaining  symbols  are  standard  (cf,  lief.  24,  pp.  173-176).  '1'he  sums  represent  the 

effect  of  discontinuities  in  the  velocity  structure  on  transmitted  (j)  and  reflected  (k)  rays. 

The  advantage  of  this  approach  is  that  accurate  results  can  be  obtained  by  evaluating  only 
three  integrals  (Eq.  ( II - 2 ) ] for  a particular  epicentral  distance  and  the  species  of  ray.  Tables 
in  Ref.  23  give  values  of  Tm(A)  for  eight  commonly  reported  body-wave  phases. 

Here,  we  purport  to  show  how  the  tabular  material  of  Ref.  23  can  be  easily  extended  to 
evaluate  ellipticity  corrections  for  surface-reflected  phases  such  as  PP,  PS,  PSS,  (ScS)^,  etc. 

Obviously,  Eq.  (II -2)  could  be  directly  applied  to  evaluation  of  ellipticity  corrections  for 
composite  phases  by  integrating  along  the  entire  ray  path,  but  this  procedure  would  be  need- 
lessly time  consuming.  For  a phase  reflected  at  a distance  Aj  and  observed  at  a distance  A^, 
tiie  integral  in  Eq.  (II -2)  can  be  represented  by  a sum  of  the  integrals 

''o  ^o  c A^ 


tiie  terms  corresponding  to  the  effect  of  the  discontinuities  are  omitted  for  brevity.  Each  of 
the  integrals  corresponds  to  the  ellipticity  corrections  accumulated  over  the  first  and  second 
legs  of  the  path,  respectively. 

Hv  changing  the  variable  O to  v i A^ . the  second  integral  now  has  the  limits  (0,  A2  - A^ ). 
With  the  use  of  some  simple  trigonometric  transformations,  it  is  easy  to  show  that 


2 

*m|,,+All=  l Amk(Al)-Xk<1') 
k-0 


(II-3) 


where 


A = 


f (1  + 3 cos  Aj) 


VI  • ,A 

~2~  sin2A1 


\fl  *2 
— sin  A1 


■s/3  • ,A 

-y  sin  2 A1 


cos  2 A. 


i sin2Ai 


■Jl  . 2 A 
~2~  Sin  A1 


- j sin  2 Aj 


1 ,,  2 A , 
j (1  f cos  A1 ) 


IJy  substituting  Eq.  (II - 3 ) into  (II -2 ),  we  find  that  the  coefficients  r corresponding  to  the  ellip 
ticity  corrections  accumulated  over  the  entire  path  of  the  ray  are 
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HI-4) 


V-V  = W - AmkIA.'VA2  V 

k=0 


r,(A,  A. ) is  the  ellipticity  correction  coefficient  for  a surface-focus,  single-leg  phase  ob 
served  at  A - A^. 

liquation  (11-4)  can  be  easily  generalized  for  an  arbitrary  number  of  surface  reflections, 
with  the  individual  legs  corresponding  to  either  1’-  or  S-wave  type. 

A.  M.  Dziewonski 


D.  CORR ESI ’ONI )ING  SHEAR  STRAIN  AND  Vp  ANOMALIES 

IN  THE  TONGA  DEEP- EARTHQUAKE  ZONE 

2 5 

In  the  previous  SATS,  Fitch  and  Jackson  presented  velocity  anomaly  patterns  computed 
from  differential  residuals  from  earthquake  locations  relative  to  a master.  To  stimulate  geo- 
physical interest,  each  pattern  and  the  corresponding  relative  locations  were  projected  onto  a 
vertical  plane  normal  to  the  strike  of  the  seismic  zone.  Each  pattern  may  be  distorted  bv  the 
defocusing  effect  of  the  high-velocity  slab,  and  each  velocity  increment  in  a given  pattern  may 
be  in  error  by,  say,  100  percent.  Of  particular  interest  is  the  pronounced  high-velocity  lobe 
in  the  down -dip  direction  of  the  deep  Tonga  seismic  zone  as  well  as  the  overall  simplicity  of 
that  pattern.  A refined  version  of  that  anomaly  pattern  is  presented  in  Fig.  II - 5.  The  master 
earthquake  is  at  a depth  of  about  6 50  km,  which  is  within  50  km  of  the  deepest  earthquakes 
located  in  this  zone.  Figure  11-6  shows  that  if  the  same  group  of  earthquakes  are  relocated 
with  respect  to  a master  chosen  from  among  the  shallower  events  in  this  group,  the  resulting 
anomaly  pattern  is  almost  perfectly  out  of  phase  with  the  pattern  for  the  deeper  master.  Fig- 
ure 11-7  shows  that  the  anomaly  in  down-dip  direction  is  negative  if  one  of  the  shallower  earth- 
quakes is  used  as  a master,  in  this  case  the  master  is  located  at  a depth  of  about  530  km.  Ap- 
parently, variations  in  near-source  velocity  structure  as  seen  from  the  shallower  master  are 
opposite  in  sense  to  those  seen  by  the  deeper  master. 

This  out-of-phase  relationship  cannot  be  the  result  of  the  different  configuration  of  the 
secondary  events  with  respect  to  each  master  or  the  velocity  increase  with  depth  expected  from 
self-compression  alone.  Both  of  these  possibilities  are  eliminated  by  computing  anomaly  pat- 
terns from  synthetic  data  that  simulate  both  the  statistics  of  the  real  data  and  the  real  distribu- 
tions of  stations  and  sources.  A Herrin  mantle  and  JIJ  core  were  chosen  for  this  simulation. 
The  patterns  for  deeper  and  shallower  masters  are  compared  in  Fig.  11-8.  The  out-of-phase 
relationship  so  obvious  in  Fig.  11-6  is  not  present.  Furthermore,  the  patterns  from  synthetic 
data  are  more  complex  than  those  from  real  data,  thus  supporting  a claim  that  the  random 
errors  in  the  real  dam  have  been  effectively  suppressed. 

The  out-of-pnase  relationship  and  other  published  evidence  strongly  suggest  that  a major 

phase  boundary  exists  somewhere  between  the  two  masters.  Jn  situ  compressional  velocities 

for  this  seismic  zone  are  high  by  5 to  10  percent  for  master  earthquakes  located  at  depths  of 

26  27 

about  600  km  and  deeper.  McGarr  showed  that  there  is  a dramatic  increase  in  the  release 
of  shear  strain  at  a depth  of  about  600  km  in  this  seismic  zone.  His  results  are  based  on  the 
conversion  of  published  values  to  seismic  moments.  The  correlation  in  depth  between  ve- 
locity and  shear-strain  anomalies  strongly  suggests  a cause-and-effect  relationship,  the 
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obvious  one  being  that  the  shear  strain  released  by  these  deep  earthquakes  is  generated  by  a 
change  to  a denser  mineral  assemblage. 

T.  .1.  Fitch 
.1.  -facksont 


E.  DETECTION  OF  MULTIPLE  EXPLOSIONS  USING  SUP  FACE  WAVES  AT  MAIO 

We  are  investigating  the  use  of  surface  waves  to  detect  multiple  explosions  from  Eastern 
Kazakh.  Although  surface  waves  are  longer  period  than  body  waves,  they  are  broader  band  in 
frequency  and  well  dispersed.  Therefore,  thev  are  well  suited  to  pulse -compression  schemes 
such  as  matched  filtering,  which  has  been  extensively  used  in  the  past1*  to  detect  interfering 
Rayleigh  waves.  At  teleseismic  distances  less  than  2 5",  the  short -period  body  waves  are  often 
quite  complex  due  to  travel-time  triplications  produced  by  the  upper-mantle  velocity  structure. 
This  often  makes  it  difficult  to  use  short-period  data  for  detecting  multiple  shots,  whereas  sur- 
face waves  at  such  distances  suffer  no  such  degradation. 

Table  11-3  lists  the  PDE  parameters  for  four  Eastern  Kazakh  explosions  with  well -recorded 
surface  waves  at  the  Meshed  SRO  (MAiO)  which  is  ~19°  from  the  epicenters.  Events  3 and  4 in 
ti  e table  have  origin  times  and  distances  which  differ  respectively  bv  5 sec  and  56  km.  Fig- 
ure II-9  shows  the  epicenters  and  great-circle  ray  paths  toward  MAIO  for  these  two  events. 


PDE 

PARAMETERS 

TABLE  11-3 

FOR  FOUR  EASTERN  KAZAKH  EXPLOSIONS 

Event 

Date 

Latitude 
(°  N) 

Longitude 

(°E) 

Origin  Time 
(GMT) 

mL 

1 

7/4/76 

49.  92 

78.95 

2:56:58.0 

5.3 

2 

8/28/76 

49.97 

79.00 

2:56:58.0 

5.8 

3 

10/29/77 

49.84 

78.  17 

3:06:59.  7 

5.5 

4 

10/29/77 

50.06 

78.91 

3:07:02.9 

5.6 

Our  scheme  for  the  detection  of  multiple  explosions  uses  least-squares  filters  to  shape  the 
one  Rayleigh  wave  to  another.  The  filter  is  interpreted  to  be  the  transfer  function  in  time  which, 
when  convolved  with  the  first  wave,  produces  the  second  wave  with  least-squares  error. 

Figure  II- 10  shows  the  Rayleigh  waveforms  (Z  component)  at  MAIO  for  Events  1 and  2 and 
the  least-squares  filter  50  points  long  which,  when  convolved  with  Event  1.  yields  Event  2 in  a 
least-squares  sense.  The  filter  was  obtained  by  solving  the  normal  equations 

R f = d (I  I - 5 ) 

where  R is  the  autocorrelation  matrix  of  Event  1,  and  3 is  the  cross-correlation  vector  between 

Events  1 and  2 evaluated  at  50  lag  values,  the  same  length  as  T.  This  equation  is  the  classic 

least-squares  Wiener-filter  equation  in  the  time  domain,  which  is  rapidly  solved  by  the  Levinson 
29 

recursion  scheme.  Because  of  the  correlations  used  in  the  equation,  this  scheme  is  analogous 
to  matched  filtering  but  more  versatile. 


t Cambridge  University,  Cambridge,  England. 


The  filter  in  Fig.  II -10  is  a veiy  simple  pulse  of  negative  amplitude  indicating  that  the  Ray- 
leigh wave  for  Event  2 is  very  similar  to  that  of  Event  1 except  for  a polarity  change,  which  is 
due  to  the  phase  difference  within  the  group  delay  envelopes  of  each  event. 

Figure  11-11  shows  the  interference  pattern  of  Rayleigh  waves  from  Events  3 and  4 com- 
pared with  Event  1.  The  least-squares  filter  to  shape  Event  1 to  the  composite  waves  of 
Events  3 and  a is  more  complex  than  the  filter  in  Fig.  II- 10,  primarily  due  to  the  two  negative 
pulses  separated  by  about  27  sec.  This  indicates  that  the  composite  Rayleigh  wave  is  well  ap- 
proximated by  a weighted  sum  of  two  Rayleigh  waves  of  Event  1. 

As  a check  on  the  stability  of  the  filter  shapes  in  Figs.  II - 1 0 and  II - 1 1 , two  extra  lags  were 
introduced  in  each  case  to  delay  the  position  of  the  negative  spikes  within  the  50-point  length 
filter.  The  new  computed  filters  are  superposed  over  the  old  ones  with  the  lags  removed,  and 
are  displayed  in  Fig.  II-i2(a-b).  The  filter  shapes  in  each  case  are  very  consistent,  and  in 
Fig.  II-  12(a)  the  two  negative  spikes  are  considered  reliable. 

The  time  delay  at  MAIO  between  Events  3 and  4 can  be  inferred  approximately  from 
Fig.  II -9 . Assuming  a group  velocity  of  ~3.0  km/sec  at  a 20-sec  period  from  Eastern  Kazakh 
to  MAIO  gives  a time  shift  of  18.  7 sec  based  on  the  difference  in  great-circle  path  lengths. 

This,  added  to  the  difference  in  origin  time  of  5.1  sec,  gives  ~24  sec  which  is  quite  close  to  the 
26-  to  27-sec  delay  between  the  negative  spikes  in  Fig.  II - 12 { a ),  considering  the  inaccuracies  in 
PDE  locations  and  origin  times. 

Unfortunately,  the  surface  waves  generated  by  this  explosion  pair  were  not  o'  rficient 
amplitude  to  be  observable  at  the  other,  more-distant  SRO  stations,  and  thus  we  are  unable  to 
unequivocally  verify  the  calculated  separation  between  surface-wave  arrivals  through  its  varia- 
tion with  azimuth.  The  short-period  recording  of  the  double  event  observed  at  MAIO  is  too 
complex  to  interpret,  and  this  may  be  a reason  for  testing  this  method  further  with  azimuthally 
distributed  surface  waves  from  multiple  earthquakes  and  explosions  when  they  occur. 

C.  W.  Frasier 
R.G.  North 


F.  SPECTRAL  CHARACTERISTICS  OF  RAYLEIGH  WAVES 

FROM  AN  AFTERSHOCK  SEQUENCE 

In  a previous  SATS,3^  we  described  an  unique  dataset  of  surface  waves  recorded  at  Meshed 
SRO  from  an  aftershock  sequence  in  the  Kuril  Islands.  The  relocated  epicenters  are  constrained 
to  an  area  of  80  x 50  km  in  extent  and,  in  addition,  over  80  percent  of  the  events  had  pP  phases 
indicating  a consistent  depth  for  the  sequence  of  45  to  55  km.  Surface-wave  magnitudes  mea- 
sured at  Meshed  range  from  4.3  to  7.3.  We  have  attempted,  without  success,  to  determine 
fault-plane  solutions  for  these  events  from  WWSSN  data:  all  we  can  say  is  that  the  few  first 
arrivals  observed  were  uniformly  compressional,  indicating  thrust  faulting.  This  is  consistent 
with  the  tectonics  of  the  area. 

We  have  considerably  refined  our  spectral-estimation  technique,  using  the  time -variable 
filtering  technique31  to  remove  both  severe  multipathing  and  the  effect  of  multiple  events.  The 
amplitude  spectra  obtained  appear  (curiously)  to  divide  into  two  distinct  types,  shown  in 
Figs.  11-13  and  11-14.  The  SRO  instrument  response  has  not  been  removed,  since  we  are  only 
interested  in  spectral  comparisons.  Type  I events  (Fig.  11-13)  have  a pronounced  spectral  hole 
at  0.025  Hz,  whereas  those  of  type  II  (Fig.  11-14)  have  two  much-less-well-developed  such  "holes." 
The  similarity  of  the  spectra  within  each  group  is  remarkable.  Because  of  the  rapid  falloff  of 
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the  SRO  instrument  response  at  periods  less  than  10  sec,  seismic  scaling  effects  will  cnly  be 
apparent  for  events  with  corner  frequencies  greater  than  0.1  Hz,  corresponding  to  those  with 
Mg  N 6.0  on  the  basis  of  the  w -square  law.  This  explains  the  almost  identical  spectra  for  the 
type  II  events,  all  of  Mg  6.3,  and  we  do  indeed  see  an  enhanced  long -period  content  for  the 
largest  type  I events. 

We  have  computed  theoretical  amplitude  spectra  for  a pure-dip  slip  event  at  50  km  depth 
in  a Gutenburg  continental  model,  and  find  that  the  spectral  hole  at  a 4j-t>ee  period  is  indeed 
predicted.  We  have  attempted  to  find  some  explanation  for  the  absence  of  this  feature  in  the 
type  II  events,  but  in  terms  of  Location,  depth,  and  (as  far  as  we  can  toll)  mechanism,  they  are 
identical  with  those  of  type  I.  We  intend  to  study  the  short -period  P-tvave  signals  from  the  two 
classes  of  events  to  see  whether  these  also  differ. 

We  have  recomputed  spectral  amplitude  ratios  fur  the  largest  events  and  find  good  agree- 
ment with  the  oj -square  (Ref.  32)  law  of  seismic  scaling.  Predicted  and  observed  ratios  are 
shown  in  Fig.  11-15.  We  also  attempted  to  measure  spectral  phase  differences  between  events, 
but  were  unsuccessful,  probably  because  the  spatial  separation  between  the  larger  events  intro- 
duces a frequency-dependent  phase  delay  for  which  we  cannot  correct  without  a more-detailed 
knowledge  of  the  dispersion  characteristics  of  the  source  region. 

R.G.  North 
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Fig.  II - 1 . Chart  showing  world  seismicity  1963-1972  and  distance  of  all  land  areas,  except  Antarctica,  from  deep-ocear. 
international  waters.  Contours  are  in  5°  intervals;  hatched  and  crosshatched  areas  show  regions  greater  than  15^  (166  5 km) 
and  20°  (2220  km),  respectively,  from  deep-ocean  areas.  Note  that  nearly  entire  earth's  surface  is  less  than  20 ° from  deep 
ocean  and  that  virtually  all  major  seismicity  zones  are  less  than  5°  from  deep  ocean. 
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Fig.  II-3.  Comparison  of  travel-time  difference  for  new  ocean-ridge  velocity  model 
and  JB  model.  Uppermost  curve  shows  differences  between  travel  times  given  in  JB 
Tables  at  an  earthquake  depth  of  33  km,  and  those  calculated  using  ray-tracing  methods 

of  Julian.2*  Small  difference  attests  to  accuracy  of  ray-tracing  method.  Four  lower 
curves  show  travel-time  differences  for  earthquakes  at  depths  (km)  indicated. 
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Fig.  II-4.  Charts  showing  relocation  of  mid-Atlantic  ridge  earthquakes;  quartered 
circles  indicate  ISC -computed  epicenters  and  their  identifying  number,  crosses  show 
relocated  positions  which  were  depth  convergent.  Open  circles  show  relocated  epi- 
center whose  depths  were  constrained  to  113  km.  Continuous  line  delineates  axes  of 
median  valley  (northerly  trending  segments)  and  transform  faults  (easterly  trending 
segments)  as  determined  from  high-resolution  bathymetric  surveys. 
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Fig.  II- 5.  Velocity  anomaly  pattern,  deep 
Tonga  master. 
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Fig.  II-6.  Comparison  o,  maly  patterns  from  deeper  and  shallower  masters. 
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Fig.  II-7.  Velocity  anomaly  pattern  shallower  Tonga  master. 


Mg.  II-8.  Comparison  of  anomaly  patterns  from  synthetic  data 
(complement  to  Fig.  II -6). 
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Fig.  II-9.  Map  of  epicenters  in  Eastern 
Kazakh  for  Events  3 and  4 in  Table  II- 3, 
taken  from  PDE  listings.  Difference  in 
great-circle  path  lengths  to  MAIO  and 
c 'igin  times  suggests  a difference  in 
group  arrival  time  at  20-sec  periods  of 
~24  sec  at  MAIO. 


Fig.II-10.  Rayleigh  waves  (Z  component) 
at  MAIO  for  Events  1 and  2 in  Table  II- 3. 
Least-squares  filter  to  shape  Event  1 to  2 
is  to  left. 


Fig.  II- 1 1.  Rayleigh  waves  (Z  component) 
at  MAIO  for  Event  1 and  sum  of  Events  3 
and  4 in  Table  II- 3.  Least-squares  filter 
to  shape  Event  1 to  combination  of  Events  3 
and  4 is  to  left. 


Fig.  11-12.  Filters  at  different  lags:  (a)  to 
shape  Event  1 to  sum  of  Events  3 and  4,  and 
(b)  to  shape  Event  1 to  2. 
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Fig.  11-13.  Rayleigh-wave  amplitude  spectra 
of  type  I,  Kuri’  Island  earthquakes  observed 
at  Meshed  SRO. 
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Fig.  11-14.  Same  as  Fig.  11-13,  type  II 
spectra. 


Fig.  11-15.  Spectral  amplitude  ratios,  observed  and  predicted,  for  two  earthquake 
pairs  of  Kuril  Island  sequence. 
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III.  COMPUTER  SYSTEMS 


A.  PDP-11  COMPUTER  FACILITY 

A steady  increase  in  computer  usage  has  put  an  increasing  strain  on  our  PDP-11  computer 
facilities.  New  hardware  purchases  and  a continuing  effort  tc  improve  existing  software  have 
been  made  to  try  to  relieve  this  strain. 

1.  New  Hardware 

The  new  hardware  acquisitions  include  a cache  memory  system,  two  new  1600-bpi  tape 
drives,  and  a 300-Mbyte  disk. 

The  cache  memory  system  is  cer^ned  to  decrease  memory  access  time  and  thereby  in- 
crease processor  throughput.  Its  installation  has  resulted  in  a 30-  to  40-percent  increase  in 
processor  throughput,  depending  on  the  type  of  load  on  the  system.  It  has  performed  flawlessly 
since  the  first  day  of  installation. 

The  new  300-Mbyte  disk  was  acquired  to  meet  the  dramatic  increase  in  the  amount  of  file 
storage  required  by  our  users.  Our  existing  disk  system  with  80  Mbytes  of  storage  was  much 
too  small  to  deal  with  the  large  seismic  data  bases  our  users  require. 

Unfortunately,  the  addition  of  the  new  disk  caused  some  unexpected  problems  in  our  overall 
system.  The  increased  load  on  the  UNIBUS  revealed  some  design  flaws  in  some  of  the  other 
peripheral  devices.  In  particular,  we  found  that  the  existing  TU10  tape  drive  controller  and  the 
VDH  ARPANET  interface  could  not  function  with  the  increased  competition  for  the  UNIBUS. 

Data  overrun  errors  on  both  devices  had  disastrous  consequences  which  made  them  almost  un- 
usable. Several  measures  have  been  taken  to  overcome  this  problem.  The  VDH  hardware  and 
software  have  been  modified  to  handle  data  overrun  errors  more  intelligently.  As  a temporary 
measure,  the  80- Mbyte  Diva  disk  system  has  been  disconnected  to  reduce  UNIBUS  contention. 
Further  improvements  in  the  software  drivers  for  the  tape  drives  and  disk  are  needed  to  com- 
pletely solve  the  pro  lem.  It  is  our  belief  that  the  newly  purchased  1600-bpi  tape  drives  will 
have  fewer  problems  because  of  improved  controller  design. 

In  addition  to  the  hardware  problems,  the  300- Mbyte  disk  has  forced  us  to  rewrite  our  file 
system  backup  procedures.  The  existing  UNIX  software  for  backing  up  file  systems  is  incapable 
of  handling  300- Mbyte  file  systems. 

The  1600-bpi  tape  drives  were  obtained  to  handle  the  increasing  tape  usage  on  our  system. 
The  planned  shutdown  of  the  PDP-7  computers  will  leave  the  PDP-11  system  as  the  only  way  to 
access  the  seismic  data  in  our  extensive  tape  library.  A number  of  hardware  malfunctions  in 
the  tape  drives  have  prevented  us  from  using  them  up  to  this  point. 

2.  Software  Work 

The  UNIX  graphics  software  has  been  modified  to  add  new  commands  and  to  speed  up  picture 
processing.  Most  of  the  changes  were  necessitated  by  the  new  seismic  data  analysis  software 
being  developed.  The  Fortran  graphics  subroutines  have  been  rewritten  in  C to  simplify  soft- 
ware maintenance.  The  translation  has  been  made  possible  by  a new  program  which  translates 
C-object  modules  to  Fortran-object  modules.  This  change  also  clears  the  way  for  a sligrt 
change  in  the  graphics  system  design.  It  will  now  be  possible  to  fully  utilize  the  inter-process 
communication  ports  without  the  limitations  put  on  them  by  our  Fortran  system.  The  graphics 
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daemons  will  now  be  created  when  a user  invokes  a graphics  command  rather  than  during  system 
startup.  The  daemons  will  also  terminate  when  the  user's  program  terminates.  In  essence,  the 
user  will  temporarily  be  the  owner  of  the  graphics  daemon  while  he  is  using  it.  There  are  sev- 
eral advantages  to  this  scheme.  There  will  be  fewer  processes  permanently  in  the  system.  It 
will  be  easier  to  specify  data  files  to  the  graphics  daemons  without  giving  complete  path  names. 
This  and  other  graphics  improvements  are  a continuing  effort. 

Heavy  usage  of  the  AH  PANIC T software  has  revealed  several  bugs  both  inside  and  outside  the 
UNIX  kernel.  A number  of  these  bugs  have  been  fixed,  but  a number  still  remain  - particularly 
two  bugs  in  the  UNIX  kernel  which  stil1  cause  occasional  systems  crashes.  In  addition,  the  FTP 
code  could  use  sor  .e  more  work  and  the  SERVER  TELNET  program  should  probably  be  com- 
pletely rewritten.  ^ gax 

B.  PDP-11  SYSTEM  SOFTWARE 

During  the  last  six  months,  we  hs.ve  worked  to  stabilize  our  system  software  and  our  main- 
tenance procedures,  and  to  prepare  the  system  for  export  to  places  such  as  the  VELA  Seisrno* 
logical  Center.  We  have  put  all  system  sources  on-line  and  carefully  checked  to  see  that  these 
sources  match  the  binaries  we  are  using.  In  a number  of  cases,  we  have  had  to  get  more  up-to- 
date  sources  because  binaries  had  been  imported  without  sources.  We  have  put  the  manual  on- 
line and  have  create  system  that  makes  it  easy  to  print  manual  sections.  However,  much 
work  needs  to  be  done  to  update  the  manual,  especially  for  system  maintenance  programs. 

We  have  also  designed  and  implemented  a set  of  maintenance  programs  called  the  Team 
Programmer  System  (TPS).  TPS  automatically  keeps  records  of  changes  made,  makes  backups 
at  the  right  moments,  and  performs  checks  when  destroying  backups.  Using  TPS,  we  expect  our 
system  maintenance  to  be  shared  with  several  other  separately  funded  projects  inside  Lincoln 
Laboratory,  thereby  distributing  the  cost.  Each  project  will  provide  one  or  more  people  who 
will  spend  typically  10  to  30  percent  of  their  time  making  system  improvements,  with  all  proj- 
ects then  sharing  these  improvements.  This  concept,  which  is  slowly  evolving,  will  be  initially 
tested  when  another  project  brings  up  a PDP-11/70,  using  our  software. 

R.  L.  Walton 
P.  A.  Neilson 


C.  DATACOMPUTER-RELATED  UNIX  SOFTWARE 

The  previous  SATS*  described  our  plans  to  develop  software  that  would  automatically  create 
Datacomputer  files  containing  associated  waveforms  for  selected  events  chosen  from  the  PESb 
(Preliminary  Event  Summary  File).  The  first  version  of  this  software  has  been  completed  and 
has  been  run  on  a test  basis  for  a period  of  two  weeks.  The  tests  revealed  that  the  software 
operated  correctly,  but  the  Datacomputer  response  was  unacceptably  slow. 

With  the  aid  of  the  user  support  services  of  the  Computer  Corporation  of  America  (CCA), 
who  are  responsible  for  maintaining  the  Datacomputer,  we  have  analyzed  the  reasons  for  this 
slowness  and  have  identified  a possible  solution.  The  first  version  of  our  software  issued  a 
separate  Datacomputer  request  for  each  waveform  segment  to  be  extracted  from  the  raw  data 
files  and  moved  to  our  own  waveform  files.  With  11  SRO  stations  presently  available,  and  a 
potential  of  1 short-period  and  3 long -period  waveforms  per  station,  this  resulted  in  up  to 
44  separate  requests  for  each  event  processed.  There  are  at  least  two  sources  of  Datacomputer 
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overhead  associated  with  each  request  - compiling  the  request  and  accessing  the  data  files  re- 
ferred to  by  each  request.  It  appears  that  the  latter  is  the  main  source  of  our  problem. 

CCA  suggested  two  ways  we  might  reduce  the  request  overhead:  (1)  to  group  multiple  re- 
quests into  a single  request,  and  (2)  to  start  using  a new  Datacomputer  feature  not  previously 
available  precompiled  requests.  Both  these  changes  involved  rewriting  the  request  in  its 
general  form  and  then  transmitting  the  parameters  of  each  waveform  transfer  to  the  request  as 
if  they  were  data.  Because  the  short-  and  long-period  raw-data  files  differ  in  format,  t he 
easiest  way  to  handle  this  is  to  group  all  the  short-period  requests  for  an  event  into  a single 
request  and  all  the  long-period  requests  for  the  event  into  a second  request.  Requests  foi  vnve- 
form  segments  that  overlap  day  boundaries  are  not  grouped  with  the  others  because  of  their 
more -complicated  format  and  the  fact  that  such  requests  do  not  occur  very  often. 

Dealing  with  these  grouped  requests  is  more  complex  than  dealing  with  single  requests. 

With  single  requests,  the  software  only  had  to  deal  with  success  or  failure  of  the  entire  request. 
With  the  group  requests,  the  software  must  be  able  to  handle  more -complicated  situations  and 
determine  which  parts  of  the  request  succeeded  and  which  parts  failed.  Recovery  after  network 
failures  is  similarly  made  more  difficult. 

When  it  became  apparent  that  major  modifications  were  needed  before  the  waveform  system 
became  operable,  we  decided  to  improve  our  interim  system  for  retrieving  SRO  data  from  the 
Datacomputer.  The  existing  program  to  retrieve  data,  DCREAD,*  requires  the  user  to  input  the 
actual  time  segment  of  data  to  be  retrieved.  In  practice,  the  seismologists  find  it  more  con- 
venient to  work  with  lists  of  events  rather  than  lists  of  stations  and  time  segments.  Therefore, 
we  have  developed  software  to  convert  event  lists  to  a properly  formatted  DCREAD  input  file. 

A program  to  do  this  for  short-period  data  (DCR1SP)  has  been  completed,  and  a long-period  ver- 
sion will  follow. 

The  event  list  used  by  DCRISP  is  in  the  format  of  the  Bulletin  lists  that  are  regularly  ob- 
tained from  the  Datacomputer  Preliminary  Event  Summary  Files.  The  second  input  to  DCRISP 
is  a file  which  contains  the  station  name  and  the  start  and  stop  times  of  each  segment  of  SRO 
short-period  data  stored  on  the  Datacomputer.  These  short-period  detection  files  are  produced 
by  a Datacomputer  program  that  we  run  each  time  new  SRO  data  are  availaole.  (Since  long- 
period  data  are  stored  continuously,  we  will  not  need  a corresponding  file  for  long-period  data.) 

Before  calling  IX'RISP,  the  user  may  edit  hij  copy  of  the  event  list  so  that  it  contains  only 
those  events  in  which  he  is  interested.  For  each  event,  DCRISP  will  compute  the  short-period 
arrival  time  at  each  SRO  station.  Before  adding  each  data  request  to  the  DCREAD  input  file,  it 
will  check  each  arrival  against  the  short-period  detection  file  and  only  retain  the  requests  for 
which  data  are  actually  available.  In  addition  to  creating  a DCREAD  input  file,  DCRISP  will 
also  pass  the  original  event  information  to  an  event  file  in  the  output  waveform  database.^  Be- 
fore actually  retrieving  the  data  from  the  Datacomputer,  the  user  may  edit  the  DCREAD  input 
file  to  further  select  out  signals  from  particular  stations  or  with  particular  distances  or  azimuths 

from  the  event.  _ , _ 

L.  J.  Turek 

P.  A.  Neilson 

D.  UNIX  SIGNAL  DISPLAY  PACKAGE 

Since  the  last  SATS*  we  have  released  the  first  version  of  the  general-purpose  interactive 
seismic  Signal  Display  Package  (SDP).  This  program  will  produce  a default  display  of  any 
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t'MX  waveform  database  and  allow  the  user  to  modify  the  display  using  a set  of  display  change 
commands. 

At  the  change  commands  described  in  the  previous  SATS  ,*  except  setting  markers  at  cursor 
locations,  plus  a few  others  not  mentioned  are  available  in  this  release  of  the  SDP.  The  change 
commands  available  in  this  version  allow  the  user  to  rescale  the  data  in  both  x and  y,  move  the 
waveforms  left,  right,  up,  or  down,  change  the  separation,  return  the  display  to  the  beginning 
of  tlie  waveform,  change  whether  the  data  should  be  displayed  clipped  or  not  clipped,  and  allows 
the  display  to  be  truncated  to  half  the  width  or  expanded  back  to  the  full  width  of  the  screen. 

The  truncate  command  is  useful  to  speed  the  plotting  process  when  only  a short  segment  of 
the  waveform  is  of  interest.  It  is  also  useful  when  using  the  special  command  which  allows  call- 
ing other  I'NIN  programs.  Truncating  the  data  to  half  the  screen  leaves  room  for  typed  input 
and  output  on  the  other  half. 

The  "post"  and  "unpost"  commands  are  also  available.  These  commands  add  pieces  of  asso- 
ciated information  to,  or  remove  them  from,  the  picture  (after  a clear  command).  The  informa- 
tion which  can  be  posted  or  unposted  includes:  the  waveform  itself,  the  date,  time,  station  name, 
instrument  code,  gram  file  name,  comment,  amplitude  scale,  timescale,  or  a zero  line.  After 
each  set  of  change  commands,  the  user  may  choose  to  have  the  changed  parts  of  the  display  over- 
written on  the  existing  display,  or  *o  have  the  screen  cleared  and  completely  replotted. 

Markers  can  also  be  displayed,  but  for  now  the  user  must  have  previously  created  his  own 
marker  file.  Work  is  still  being  done  to  allow  cursors  to  be  displayed,  and  markers  to  be  placed 
at  the  cursors.  This  feature  will  be  included  when  the  second  version  of  the  SDP  is  released. 

Two  special  commands  have  been  included  in  this  version.  The  first,  "goback",  returns  the 
display  to  the  previous  picture.  Only  one  previous  display  is  saved.  The  second,  "1",  allows 
the  user  to  call  any  of  his  own  programs  or  UNIX  system  programs  to  do  waveform  processing, 
make  nonstandard  modifications  to  the  display,  or  copy  and  save  current  display  command  files 
and  then  return  to  the  display  program. 

By  releasing  the  current  version  in  its  incomplete  state,  we  have  been  able  to  get  some  feed- 
back from  the  seismologists  for  improvements  and  programs  to  be  added  to  the  display  package. 

When  the  cursor  opt1* on  is  available,  we  will  be  able  to  add  features  that  utilize  it,  such  as 
a command  to  align  the  waveforms  according  to  the  cursor  and  programs  which  work  on  window's 
of  data  which  have  been  defined  by  the  cursor. 

Some  special  two-pass  commands  are  planned  which  would  delete  or  add  waveforms  to  the 
display,  reorder  the  existing  waveforms  (they  are  currently  displayed  in  the  same  order  in  which 
they  appear  in  the  gram  index  file^),  or  make  a display  based  only  on  "posted"  waveforms. 

One  of  the  new  programs  planned  is  a 3-component  rotate  program  which  would  locate  grams 
that  are  3 components  of  the  same  instrument  and  convert  them  to  vertical,  radial,  and  trans- 
verse components.  Another  program  will  locate  grams  that  are  3 components  of  the  same  instru 
ment  and  assign  the  same  vertical  gain  to  each  of  them. 

We  are  also  trying  to  make  general-purpose  programs  from  programs  previously  written  by 
the  seismologists  for  their  own  specific  projects.  Among  these  are  a filtering  program,  a dec- 
imation program,  and  an  equal  gain  program  to  set  the  gains  of  all  waveforms  to  an  equal  value. 

To  help  us  in  writing  these  general-purpose  programs,  we  are  compiling  a package  of 
C-language  and  Fortran  subroutines  to  make  working  with  the  database  files  easier.  We  have  al- 
ready provided  a set  of  UNIX  commandr  that  can  manipulate  all  files  in  a waveform  database  with 


70 


a single  command.  WRM  will  delete  a waveform  database,  WMV  will  rename  a database  (which 
involves  modifying  every  file  that  contains  references  to  the  database  name),  and  WDIJP  will 
duplicate  a database  under  a new  name,  with  the  option  of  forming  links  to  ihe  gram  files  in  order 
to  save  disk  space.  M ,,, 


L.  .1.  Tur-k 
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GLOSSARY 


ARPANET 

DARPA  Computer  Network 

CCA 

CTBT 

Computer  Corporation  of  America 
Comprehensive  Nuclear  Test  Ban  Treaty 

FTP 

File  Transfer  Protocol 

ISC 

International  Seismological  Center 

■IB 

Jeffreys -Bullen  (Tables) 

LASA 

LRSM 

Large  Aperture  Seismic  Array  (Montana) 
Long  Range  Seismic  Measurements 

NORSAR 

NTS 

Norwegian  Seismic  Array 

Nevada  Test  Site 

OBS 

Ocean-Bottom  Seismograph 

PDE 

PESF 

PNE 

Preliminary  Determination  of  Epicenters 
Preliminary  Event  Summary  File 

Peaceful  Nuclear  Explosion 

SATS 

SDP 

SNR 

SRO 

Semiannual  Technical  Summary 

Signal  Display  Package 

Signal-to-Noise  Ratio 

Seismic  Research  Observatory 

TPS 

Team  Programmer  System 

USGS 

U.S.  Geological  Survey 

VDH 

Very  Distant  Host  Interface 

WWSSN 

World-Wide  Standard  Seismograph  Network 
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